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Abstract 


Doppler shift measurements of the Na Dj absorption line have revealed solar oscilla- 
tions in a new regime of frequency and wavenumber. Oscillations of vertical velocities 
in the temperature minimum and low chromosphere of the Sun are observed with fre- 
quencies ranging up to 9.5 mHz. There is no evidence for chromospheric modes of 
3 minute period. This indicates that the chromosphere does not form a good cavity 
for acoustic waves. 

The fundamental-modes appear with wavenumbers up to 5.57 Mm -1 (equivalent 
spherical harmonic degree, 3877). The frequencies lie below the predicted values at 
wavenumbers above 1 Mm” 1 . The values are in agreement with previous measure- 
ments that exist for wavenumbers up to 2.67 Mm -1 . 

Spatial maps of velocity power show that high wavenumber oscillations are sup- 
pressed in active regions. The shape of the power depression indicates that wave 
motion is affected in the layer of atmosphere where the measurement is made. The 
/-modes are suppressed in the same way as p-modes, indicating that the mechanism 
for wave suppression affects velocity fluctuations. Mode frequencies are not affected 
by the magnetic fields by more than 50 /tHz, the precision of the measurement. 
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Chapter 1 
Introduction 


Solar oscillations were first observed by Leighton (1960), and were explained by Ulrich 
(1970) and Leibacher & Stein (1971) as being the superposition of many global stand- 
ing waves. They have subsequently become an important tool for the determination 
solar structure. Measurements from helioseismology include the depth of the con- 
vection zone (Gough 1976), helium abundance in the convection zone (Christensen- 
Dalsgaard & Gough 1980), and internal rotation rate (Duvall et al. 1984; and others). 

In the past, observational limitations have restricted helioseismology to long wave- 
length oscillations (~ 4 Mm). New instruments and techniques are beginning to make 
shorter wavelengths accessible to study. This opens the door to the exploration of 
structure at much smaller spatial scales. Stratification of the atmosphere, and vertical 
structure of magnetic regions are becoming active areas of investigation. 

The term solar oscillations conventionally refers to two types of waves. Because 
they exist in a resonant cavity, these waves are also called modes (eigenmodes of 
the cavity). Compressional or acoustic waves (p-modes) are longitudinal waves with 
pressure as the restoring force. They have at least one nodal surface in the radial 
direction. Gravity waxes have gravity as the restoring force. A special case of these 
are the surface waves. Because they have no radial nodes, surface waves are called 
fundamental modes (/-modes). These transverse waves propagate at the interface 
between two media. The remainder of the gravity waves are called p-modes. To date, 
p-modes have not been unambiguously detected, and they will not be considered here. 
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1.1 Wave Properties 


1.1.1 Acoustic Waves 

Sound waves can propagate vertically at frequencies greater than the acoustic cutoff 
or critical cutoff frequency u> c . At lower frequency, a disturbance moves the whole 
atmosphere up and down in phase. Hence, waves do not propagate. At low horizontal 
wavenumber k ( k <C ^tj). the acoustic cutoff can be expressed in terms of the adia- 
batic sound speed c and 1 he density scale height H . For an isothermal atmosphere, 
the scale height is given l>v 

c 2 

H = ~ (1.1) 

19 ’ 

and, the acoustic cutoff, by 


U) c = - 


2 H 


( 1 . 2 ) 


where g is the acceleration due to gravity; and, 7 is the adiabatic index (Lamb 
1932). Taking c = 7.1 km/s, g — 27S m/s 2 , and 7 = | in the photosphere, we find 
H = 110 km and u> c /2tt = 5.2 mHz for k <c 4.5 Mm -1 . 

In geneial, the acoustic cutoff depends on wavenumber. The relationship is given 
by (Kumar & Goldreich 19S9) 


ut-' 


= (oT/) 2 T It 1 + (^tf) 2 ] + y[l + (2^p_g(2^) 2 


(1.3) 


Figure 1.1 shows <jj c as a function of k. 

Oscillations with wavelength greater than a few Mm propagate around the sun 
without excessive dissipation. A periodic boundary condition is, therefore, imposed 
on these waves. The wave equation is separable into radial and angular parts and has 
eigensolutions of the form, 

$ nlm (r, 0 , 0 ,t) = R nl (r) >7"(M) e-^ (1.4) 

n = 1.2,3 / = 0, 1,2 m = -/,-/+ 1 ,...,/ 

Short wavelength oscillations are dissipated before traveling around the sun. More- 
over, at their small spatial scales, the medium can be approximated as plane-parallel. 
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Plane waves, therefore, provide a more natural basis set for the horizontal eigenfunc- 
tions. 

*»i'(r,x,y,t) = R nk (r)e i ' i - i -“ t ' (1.5) 

The horizontal wavenumber h = |A-| can be related to the spherical harmonic degree 
/ as follows. 

,//(/ + 1 ) 

k = ( 1 - 6 ) 

“sun 

Ray paths for a p-mode are diagrammed in figure 1.2. As waves travel downward, 
they encounter higher temperatures and an increasing sound speed. This causes them 
to refract back upward where they encounter increasing values of u> c . When the wave 
can no longer propagate < *^ c ), it is reflected back down. The resonance condition 
is met when crests of downgoing waves are superposed on crests of refracted waves 
returning to the surface. For waves trapped in an adiabatic and polytropic layer 
(polytropic index = 4), this condition gives a dispersion relation of approximate 
form, 
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Figure 1.2: Propagation of an acoustic wave in a plane-parallel atmosphere. The wave 
is refracted as it propagates downward, and reflected at the surface. The maximum 
depth reached by the wave is called the lower turning point. From Gough (1983). 


w ' 2 = (£ + 0 (jk n = 1,2,3,... (1.7) 

where g is the acceleration due to gravity at the surface. 

This dispersion relation can be measured directly from the ridges of power in a 
spectral decomposition (a plot ot power vs. u.' and k) of solar velocities or intensities. 
The exact relation depends on the variation of sound speed with depth. 

The lower turning point d ot the waves is approximately given by, 


d = 


3(n + 2) 
4k 


( 1 . 8 ) 


The excitation source of the p-modes is thought to be turbulent convection at 
the top of the convection zone. The efficiency of acoustic emission from convection 
is proportional to approximately the eighth power of the Mach number (Goldreich 
& Kumar 1988). Hence, only a few of the granules on the high velocity tail of the 
distribution are expected to contribute to the acoustic noise driving the p-modes 
(Brown 1991 ). 


1.1.2 Surface Waves 

Surface waves can propagate on the interface between convectively stable media (e.g., 
air over water). On the Sun, the transition region between chromosphere and corona 
provides such a boundary (.Jones 1969). The model of Avrett (1992) describes the 
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Figure 1.3: Temperature and density of the solar atmosphere. From Avrett (1992). 

temperature and density discontinuities shown in figure 1.3. For purely vertical den- 
sity stratification of a compressible fluid, surface waves obey the dispersion relation 
(Whitney 1958), 

u > 2 =gk (1.9) 

where g is the local acceleration due to gravity. 

The simple form of this relation has prompted the remark that, “There is no 
interesting physics in the /- modes.' 1 It wasn’t until /-modes were found not to obey the 
above relation at high wavenumber that this conclusion was reexamined. Libbrecht, 
Woodard &; Kaufman (1990) found /-mode frequencies to lie below the predicted 
values for wavenumbers above 1.4 Mm" 1 . At least three possible explanations exist, 
though numerical predictions have not yet been offered. 

Campbell & Roberts (1989) and Evans h Roberts (1990, 1991) predict that a 
horizontal magnetic field in the chromosphere would increase the frequencies of the 
/-modes by increasing the restoring force on the surface waves. The sign of the 
frequency shift may reverse at high wavenumber, however, in a more realistic model 
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of the chromosphere with an appropriate magnetic field structure (Roberts 1991). 

Rosenthal (1992) finds that evanescent waves dissipated in the corona would cause 
power to leak out ol the /-modes. This would lower the frequency of the modes at a 
given wavenumber by an amount dependent on the temperature difference between 
chromosphere and corona. 

Finally, structures in the atmosphere, such as flux tubes and convection cells, 
violate the assumption of horizontal homogeneity. It is not known whether this should 
increase or decrease the mode frequencies (Gough 1992). 


1.2 High Frequency Modes 

For the purpose oi discussion here, t he high frequency regime includes modes with 
u; > u? c , where u> c /2x = 5.2 mHz. This value of u> c is the photospheric value for low 
wavenumbers. Because waves with ^ can propagate, they are not trapped by 

the photosphere. So, in the model illustrated by figure 1.2, there should be no upper 
reflecting surface at frequencies > 5.2 mHz. Removal of the top surface of the cavity 
eliminates the possibility of resonance, and hence, the mode structure. 

It was a surprise, then ‘lore, when the p-mode ridges were found to continue above 
the acoustic cutoff (Jeffries el al. 1988; Libbrecht 1988; Duvall et al 1991). Two expla- 
nations have been put forward lor the origin of the high frequency ridges. Balmforth 
& Gough (1990) propose that strong reflection from the base of the corona defines a 
chromospheric cavity which couples to the subphotospheric cavity. Waves with fre- 
quency high enough to pass through the photosphere would still be trapped by the 
corona. 

The presence of a chromospheric cavity has further consequences, however. The 
chromosphere is of order one wavelength in thickness for 5-minute modes; and waves 
of appropriate wavelength A,, will resonate in this cavity. The p-modes propagate 
almost vertically in the atmosphere (angle of incidence ^ 0). Thus, A r — cju r , 
independent of horizontal wavenumber. The chromospheric modes should, therefore, 
have a dispersion relation, 


yjj = const. 


( 1 . 10 ) 
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Figure 1.4: Mode frequencies vs. wavenumber with a chromospheric cavity present. 
From Ulrich k Rhodes (1977). 

Ulrich k Rhodes (1977) found that the p-mode ridges would merge with a set of 
chromospheric modes at about 5.5 mHz. See figure 1.4. The p-modes would continue 
on above the 5.5 mHz band with an “avoided crossing” between the ridges above and 
below the chromospheric resonance. 

Dame, Gouttebroze k Malherbe (1984) claim to have seen the chromospheric 
modes in Ca II K-line observations. Their power spectra are very noisy, however, and 
it is difficult to discern any definite ridges of power. Deming et al. (1986) find a peak 
in the one-dimensional power spectrum (power versus u;) at 4.3 mHz. Its significance 
is difficult to assess, though, without knowledge of the wavenumber dependence of the 
power. Woodard k Libbrecht (1991) do not find this peak in similar one-dimensional 
power spectra. 

The other explanation for the high frequency ridges was offered by Kumar et al. 
(1990). They propose that waves are emitted isotropically from a source near the top 
of the convection zone. Waves traveling directly toward an overlying observed surface 
interfere (in Fourier space) with waves refracted upward from below to produce “mock 
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modes. ’ It must be stressed that this interference is not physical. In fact, the directly 
emerging wave trains and refracted wave trains may arrive at the surface at different 
times and in different places. The basis functions, used to represent the waves, can 
only distinguish time, modulo 2ir/u>, and position, modulo 27 t/£. Frequencies of the 
mock modes can be used t o determine the depth of the acoustic source layer (Kumar 
k Lu 1991). 

The two models for high frequency waves can be distinguished by observations 
that resolve the ridge structure in the 5.5 mHz frequency range. 

1.3 High Wavenumber Modes 

Most of the literature in helioseismology refers to the (spherical harmonic) degree 
/ of a mode rather than the wavenumber k. By convention, low degree refers to 
0 < / < 4, intermediate degree, to 4 < / < 100, and high degree to 100 < /. The 
oscillations presented in this dissertation have 500 < / < 3877. By analogy with 
the electromagnetic spectrum, they should, perhaps, be called very high or ultrahigh 
degree, but I will stick to the term, high wavenumber. 

Previous observations oi the 7-modes exist up to wavenumbers of 2.67 Mm -1 
(Libbrecht, Woodard Sz Kaufman 1990). As mentioned above, the frequencies do 
not agree with equation 1.9 at wavenumbers above 1.4 Mm -1 . The variation of g 
with height accounts (or only a 0.5% modification to equation 1.9. This is an order 
of magnitude less than the observed effect. It would be desirable to confirm the 
observational result and determine the trend at higher wavenumbers. 

The size of structures that can be resolved by oscillations is roughly equal to 
the wavelength of the modes. For active regions, spatial scales range from 200 km 
(flux tubes) to 25 Mm (large sunspots). The corresponding wavenumbers range from 
0.25 to 30 Mm 1 , Although the high end of this range is not attainable with any 
current instrument, there is still much room for investigation. 

Woods Sz Cram (1981) have reported a suppression of oscillation amplitudes in 
plage as compared to quiet re gions. Tarbell et al. (1988) and Title et al. (1992) find a 
suppression of both convective motions and wave amplitudes. The suppression of the 
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oscillations coincides with the boundary of the photospheric magnetic field to within 
one wavelength. 

Braun, Duvall & LaBonte (1987, 1988) and Braun et ai (1992b) have reported 
absorption of p-mocle power by sunspots. The mechanism for this is not understood. 
Resonant absorption (Ilollweg 1988; Chitre & Davila 1991), scattering to short wave- 
length (Bogdan Sz Zweibel 1987), and enhanced viscous dissipation (LaBonte h Ryu- 
tova 1992) are among the explanations offered. It is by no means clear whether 
only propagating acoustic waves below the photosphere are absorbed, or whether 
evanescent waves are also affected. A spatial map of oscillation amplitudes above 
the photosphere, where the magnetic fields are different, would help untangle the two 
effects. 

Global oscillations have been found to have frequency shifts correlated with solar 
activity (Palle, Regulo &z Roca Cortes 1989; Woodard et ai 1991). This may be due 
to a change in the shape of the cavity surface with sunspots present. Waves of a 
given frequency may have different wavenumbers within magnetic and nonmagnetic 
regions. Global measurements would then find a weighted average of the two disper- 
sion relations. Comparison of high resolution observations in quiet sun and active 
regions could resolve this issue by direct measurement of dispersion relations in the 
two types of regions. 

1.4 Height of Observations 

Oscillations can be observed as fluctuations in intensity and line-of-sight velocity. 
Both p-modes and /-modes produce motions of gas particles, and cause Doppler 
shifts of absorption lines. Pressure fluctuations due to p-modes produce tempera- 
ture/intensity fluctuations. The /-modes also produce intensity fluctuations, though 
the reason is not well understood. It may be that the ~ 100 km vertical displacements 
of parcels of gas affect the process of absorption line formation. 

Different parts of an absorption line are formed in different layers of the atmo- 
sphere. Hence, the height at which the velocities and intensities are determined 
depends on the positions in the absorption line at which measurements are made. In 
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Figure 1.5: Na D line profiles on left hand scale. Contribution functions for emergent 
intensity on right hand scale. Thick lines show the maximum of the contribution 
function at each wavelengt h. Thin lines show 1/3 of the maximum of the contribution 
function. From Schleicher (1976). 

the core of deep lines such as Na D lf very little continuum light emerges from the 
photosphere. The height of a measurement, therefore, depends in the contribution 
function for the emergent line emission (Gurtovenko, Sheminova & Sarychev 1991; 
Schleicher 1992). This function is shown in figure 1.5. 

The original intention of the observations described in chapter 4 was to make one 
set of measurements in the core and one in the wings of the Na Di line. Problems 
described in section 4.1 resulted in a measurement weighted mostly in the core (600- 
850 km), but with some contribution from below the temperature minimum ( 1 50— 
500 km). 



Chapter 2 

Measurement Techniques 


Producing power spectra from uncalibrated images involves several stages of data 
analysis and reduction. This chapter reviews some general procedures, standard and 
nonstandard, for measuring line-of-sight. components of velocity and magnetic field, 
and for determining wave amplitudes and phases. The details of calibration to remove 
artifacts of the instrument and atmosphere are left to chapter 5. 

The line-of-sight velocity is found trom the Doppler shift of an absorption line. 
The Doppler shift, in turn, is determined by measuring intensity as a function of 
wavelength at a series of fixed- wavelength positions in the absorption line. The more 
points in the measurement, the better the velocity estimate will be. Unfortunately, 
more points require proportionately more time to make the measurements, and geo- 
metrically more time to analyze them. 

In this chapter are described two new techniques for making velocity images (dopp- 
lergrams) and line-center intensity measurements using three-point and four-point 
methods. The four-point method was developed for the Solar Oscillations Investi- 
gation / Michelson Doppler Image r (SOI/MDI) experiment for the Solar and Helio- 
spheric Observatory (SOI 10) spacecraft. It was not applied to the data used for this 
dissertation. 

Images of line-of-sight magnetic field strength (magnetograms) are an important 
tool in this study. The tried and tested "two-point magnetogram’' method is described 
below. 


11 
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Figure 2.1: Measurements of the Na Di line 

Finally, the Fourier transform is essential to the study of waves propagating in 
a linear medium. Procedures used to extract information about the oscillations are 
described in the last section of this chapter. 


2.1 Three-Point Measurements of Line Parame- 
ters 

The three-point method for finding Doppler shifts and line center intensities assumes 
a parabolic absorption line profile. This assumption is valid when the three intensity 
measurements are positioned relatively close to the line core. Figure 2.1 shows the 
measurements schematically for the Na D! (5896 A) line. The Doppler shift and line 
center intensity are the j and y coordinates of the vertex, respectively. 
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2.1.1 Algebraic Solution 

Assuming a parabolic lino profile 


/(A) = / 0 + 


A — Aq 

const. 


( 2 . 1 ) 


we can solve for the vertex (Ao« /o) m terms of three intensities measured at 

wavelengths Ai,A 2 ,A 3 , respectively. 

Ai ci + A3C3 


Au — 


h) — I2 + 


C\ + C 3 

(Vi - 6 3 c 3 ) 


2 


1^1 ^3 (^'3 “ ^l)( c 3 — c l) 

where A! = \(X\ + A 2 ), A 3 = ^(A 3 + A 2 ), 

< s 1 = A i — A 2 . S 3 — A 3 — A 2 , 

< i = 6 \{I 3 - / 2 ), and c 3 = <5 3 (/i - h). 


( 2 . 2 ) 

(2.3) 


2.1.2 Error Estimates 

The errors on the measured quantities /i,/ 2 ,/ 3 are due to limitations of the instru- 
ment. The transmission ol the tunable filter varies with wavelength as shown in 
figure 3.8. The phase of each Fourier component ol this modulation is independent, 
and drifts with time. I 11 addition, the wavelengths tuned by the filter are quantized 
in units of the smallest step size (3.14 mA). These sources of error were used in a 
simulation to estimate errors tor the velocity and intensity. 

The simulation assumed a parabolic line profile with parameters measured from 
the data. A range ol velocities ol ±2000 m/s with a flat distribution, and a range of 
line-center intensities of ±8% with a Gaussian distribution were used in the model. 
Intensity “measurements I\, / 2 , / 3 were modeled assuming a random jitter of up to 
±1.57 mA (flat distribution) in the location of each measurement, and assuming 
phases for the transmission modulation which varied randomly between sets of 3 
filt erg rams, but not within a set. The velocities and intensities were then calculated 
using equations 2.2 and 2.3. I he errors had approximately Gaussian distributions 
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with standard deviations. 


n,. — 300 m/s, and 07 = 7.7%. 


(2.4) 


2.2 Four-Point Measurements of Line Parame- 
ters 


An algorithm was developed for the SOI/MDI experiment to determine Doppler shifts 
using four measured intensities. The measurements are used to solve a system of si- 
multaneous equations for the four absorption line parameters. Since computing time 
must be kept to an absolute minimum, the nonlinear equations are solved (numeri- 
cally) using a look-up table approach. 

The model included only the Lyot and Michelson filter elements, and the following 
parameters were used: 


• Lyot filter: profile Gaussian 

fwhm 500 mA (see section 2.2.3) 

• Michelson elements (2): profile — - cos’ 2 (A; A) cos 2 (2 A A) 

fwhm — 105 mA 


where k 


Ju'hr, 


COS 


-1 2-0 




• Ni I (6768 A) L ine: profile — Gaussian 

fwhm — 130 mA 
depth --- 58% 


• Measurement Spacing — 85 mA (symmetric about nominal line center) 


See figure 2.2 for illustration. 

Four measurements are made at increasingly longer wavelengths across the ab- 
sorption line, to yield intensities l\ through / 4 . The ratios, a and /?, are formed as 
shown: 


n 


h~h 

U-h 


2.5 
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Figure 2.2: Spectral profiles. 


3 = 


1 7 ~ /, 


(2.6) 


/ 1-/4 

The quantities a and 3 are independent of CCD dark current and gain variations. 
They are also independent of continuum intensity and line-depth if a suitable choice 
is made for the Lyot filter profile (see section 2.2.3). 

We can solve exactly for the Doppler shift Ao, line depth d, and continuum inten- 
sity Iq if we assume a cosine shaped absorption line with measurement spacing equal 
to the hwhm: 

/(A) = lo 

Solving, we find, 


wt 


d ( tt(A - Ao) ^ 

» [* - ( C0S fwkm ‘jj 

(2.7) 

II 

£ 

3 

1 

n 

O 

c+ 

1 

(2.8) 

/o = I (/ 1 + / 3 -s) 

(2.9) 

d = T 

(2.10) 

I 0 

s = yjih - h ) 2 + (h - I2) 2 
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Figure 2.3: Ratio a vs. velocity for width increments: 0, 10, 20, 30%. 


Note that a has a one-to-one correspondence with the ratio r = 
which appears in theSOI/MDI proposal (i.e., r = (a — l)/(a+l)). Hence determining 
the velocity from a alone is equivalent to determining it from tan -1 
For a Gaussian line-profile, 


r. 


/(A) = / 0 


1 — de 


iff* 


( 2 . 11 ) 


where a = fwhm/ 2.36 

The ratios a and (3 are plotted as functions of velocity for different line-widths in 
figures 2.3 and 2.4, respectively. It is apparent that the velocity determined from 
a alone would be in error if the calibration table were not made for the “right” 
line- width. 

This caveat can be removed by determining the velocity v from both a and /3: 
v = v(a,{3). An algorithm for doing this is outlined below. 
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Figure 2.4: Ratio 3 vs. velocity for width increments: 0, 10, 20, 30%. 

2.2.1 Implementing the Algorithm 

Look-Up Table Approach 

It is possible in principle to compute velocity by using a two-dimensional look-up table 
in a and (3 (i.e., v — t>(o,/?)). This, however, would require an inordinate amount 
of memory ((12 bits in a) x (9 bits in 0) = 2M words) to achieve an accuracy of a 
few m/s. A better strategy is to take advantage of the fact that a, alone, provides a 
pretty good estimate of velocity; and write: 

v{a,0) = vi{a) + v 2 {a,0) (2.12) 

where V 2 is a correction to the velocity i>j and, hence, does not vary with a as quickly 
as does v\. 

This approach can achieve an accuracy similar to the previous approach but with 
much less memory overhead. The table, Uj(a), requires a 12 bit address; and the 
table, vi(a,0), requires 7 bits in a by 9 bits in 0 (this has not been optimized); for 
a total of 66k words. 
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Domain and Range of o and ;i 

The ratio a has a pole in the negative velocity-domain (i.e., a red-shifted line). To 
avoid this, a test is made to determine the sign of the velocity. For positive velocities: 
a = • For negative velocities: a = h ~ h 

The velocity is negative if and only if: I x — I 3 > I 4 — I 2 (2.13) 


As can be seen from figure 2.3, a covers a more-or-less symmetric range about 
zero. Therefore, it can be scaled-up to provide an address for a look-up table without 
having to subtract an offset. 

The function T is symmetric in the velocity domain and is continuous for sym- 
metric line profiles. For asymmet nc lines (and noisy lines) 3 has a pole near v = 0. 
In practice this does not cause problems because, near v — 0, v is independent of (3 
(i.e., v 2 (a,3) % 0 for small r). An offset may have to be subtracted from j3 when 
scaling it for use as an address. 

Finally, note that: 


h-l 3 _ (A - A) - (A - A) 
A - A (A - A) 


(2.14) 


Hence, both a and 3 can be obtained from the differences (/j — / 3 ), ( I x — / 4 ), and 

(A - A). 


The Algorithm 


It is assumed here that the differences: A x = I x — / 3 , ,4 2 = A — A, and A 3 = / 4 — I 2 
have been accumulated in image buffers. 

For all pixels: 


Calculate 3 and sign of velocity: 

numerator — A { — , 4 3 

{ — 1 , if numerator > 0; 

+ 1. otherwise. 

Q numt nilnr ^ » 

■J - y- ' i 
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Calculate a: 



Calculate velocity: 


it sign > 0: 
otherwise. 


u = sign ■ [ 1 1 ( C 2 o ) + V 2(^3 o , ( \ S)} 


are constants: \\ and V 2 are look-up tables. 

It should be noted that noise in the data, may cause a “division by zero" error in 
the algorithm. 


Constructing the Calibration Tables 

The following procedures are used to construct the calibration tables Vj and V 2 . 
Parameters used for this study appear in square brackets. 

The table \ \ ( o ) is made tor a line-profile of nominal width [130 mA] as follows: 

1. Calculate a for a representative number [44] ot discreet (evenly spaced) 
velocities ranging between some number slightly less than zero [-40 m/s] 
and some number slightly larger than the largest velocity expected from 
the data [4500 m/s]. 

2. Invert the table of o versus v to get a table of v versus a. 

3. Interpolate (cubic spline) between the endpoints to get Vi (a), a. table of 
velocities for an evenly spaced abscissa a. 

The table V- 2 (o\/3) is made for the same range of velocities as is V x \ and, for a 
given range oi widths [100 to 205 mA]. The minimum width should not be less the 
measurement spacing [85 mA]. 

There are two "unphysicar regions in the table V 2 * The highest values of a cannot 
occur with the highest values ot 3. and the lowest values of ot cannot occur with the 
lowest values of 3 . These portions of the table will not normally be addressed by the 
Doppler algorithm. 

V 2 is calculated as follows: 
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1. Calculate o and :i for a representative grid of velocities [22] and widths 

[14]- 

2. Interpolate (spline) to get a much finer grid [512 in velocity by 256 in 
width], 

3. For each "real" velocity v [512 of them] calculate the velocity correction 

t> 2 : 

t '2 = v -V^a) 

using the table \ ] ( o I computed above. 

1. Sort theorderc'd triplets (r 2 . a, 3) into a 2-dimensional table V 2 ( o , 3), keep- 
ing track of how many times each bin gets filled. 

5. Fix uncooperative Iritis: 

a) If a bin accumulated more than one value of v 2 , normalize to get 
the "average” c> for that ( a . 3 ) . 

b) If a bin in a "physical” region did not get filled [there were 6 of 
them], interpolate [linearly] using adjacent bins. 

c) Fill "uuphysicar bins with the value of v 2 from the nearest “physi- 
cal bin. 1 his is done because noise in a measurement may produce 
an ordered pair (o,T) which does not lie quite where it should in 
the table. 

Tables Vj and V 2 were made in units of m/s, and rounded to the nearest integer. 

2.2.2 Spectral Simulation 

The instrument and Doppler algorithm were simulated in order to determine the 
magnitude of statistical and systematic errors. The simulation program stored the 
line and transmission profiles in single-precision 512-element tables covering the range 
±800 mA about nominal line-center. The intensities, / 1 , . . . , I 4 , were calculated for a 
discrete set of velocities as follows: 


1. Offset the index of the line-profile table to simulate a Doppler shift. 
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2. Take the product of the transmission-profile and the offset line-profile. 

3. Integrate [sum] over wavelength to get "measured" intensity. 

The Rienian n interval and limits of integration were varied to verify that they 
did not affect the results. 

The error estimates are described in the next two sections. 


2.2.3 Sources of Systematic Error 

Line-Depth Variations 

The ratio a will be sensitive to variations in line depth if measurements I\ and 1$ (or 
/ 4 and / 2 ) do not integrate the same amount ot continuum. 

Figure 2.5 shows the transmission through the Lyot and Michelson filters as a 
function of wavelength for the first and third measurements. The ratio of integrals of 
these curves is plotted as a function of Lyot filter fwhm (assuming a Gaussian Lyot- 
profile) in figure 2.6. The two integrals are equal when the Lyot fwhm is ~ 625 mA. 
They are different by 2V( when the Jwlun is 500 mA. 

This 2% difference causes a velocity error t v which depends on line- depth as shown 
in figure 2.7. If the Lyot fwhm were 400 mA. t v would be larger by a factor of 4. 

Lyot Filter Drift 

If the Lvot filter profile drifts in wavelength space, the four measured intensities may 
not integrate the same amount ot continuum. Ihis would give erroneous values for ot 
and f3 and, hence, v. To eliminate this source ol error the transmission profile ti yo t(\) 
should satisfy: 

J */,*/( A - A\) /'„„>/,( A)rfA = constant i = 1 4 (2.15) 

for some small range d\ over which the Lyot filter might drift; t l mich (X) is the tians- 
mission of the Michelson filter at wavelength setting i ; 

For a Gaussian Lvot-profile, it turns out, this condition is more or less satisfied 
for the same fwhm that minimizes e v due to line-depth variations. The reason can be 
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Figure 2.5: Transmission product for Lyot and Michelson filters vs. wavelength. 1 1 
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Figure 2.6: Transmission ratio f 1 3 d\f f t\ dX as a function of fwhm of the Lyot filter. 
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Figure 2.7: Systematic error e„ vs. velocity for line-depth increments 

- 20 ,- 10 , 0 , 10 , 20 %. 


seen from figure 2.5. For transmission profile, if the Lyot profile shifts to the left 
(right), any increase (decrease) in the area under the main-lobe is compensated by a 
decrease (increase) in the area under the largest side-lobe. Transmission profile t$ is 
unchanged to first order since it lies close to the maximum of the Lyot- profile. 
Figure 2.8 shows e„ for 8 mA offsets of the Lyot filter. 

Line-Width Variations 

As stated above, if a alone is used to compute v, an error will be introduced which 
depends on the line-width. Figure 2.9 shows the error incurred if the line-width is 10, 
20, and 30% larger than the nominal width for which the calibration table was made. 
If both oc and 3 are used, the error is reduced to a few m/s as shown in figure 2.10. 
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Figure 2.8: Systematic error t v vs. velocity for Lyot drifts: —8,0.8 mA. 
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Figure 2.9: Systematic error t L , vs. velocity for line-width increments: 0, 10, 20, 30% 
(a-only method). 
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Figure 2.10: Systematic error e u vs. velocity for line-width increments: 0, 10, 20, 30% 
(a, j3 method). 

Line- Asymmetries 

The effect of line-asymmetry was explored by using a line-profile with a Gaussian- 
shaped bisector as shown in figure 2.11. It was found that the depth at which the 
velocity is measured depends on the velocity, but is independent of bisector shape. 
The velocity dependence is plotted for different bisector shapes in figure 2.12 (a-only 
method), and figure 2.12 (a. 3 method). Note that e v does depend on the precise 
shape of the bisector, but its magnitude cannot be determined until a measurement 
is made of the bisector shape at the spatial resolution of the MDI instrument. 

2.2.4 Statistical Error 

Statistical errors were found by a standard propagation of error calculation. Deriva- 
tives of the functions ui(cv) and v? (a, 3) were approximated by finite differences. 
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velocity (m/s) 

Figure 2.13: Depth of measured velocity vs. velocity for bisector shape hwhm: 

—7.5, . . . , 7.5 mA (a. 3 method) 


For a sinusoidal line-profile with delta-function transmission- profiles, the statisti- 
cal error a v can be written down in a form like that shown in the SOI/MDI proposal. 


2cs 1 y/2 -d 

7rA v/7o d 


(2.16) 


where, s is the measurement spacing; d, the line depth; and, 7o, the number of photons 
collected at continuum intensity (half the well depth for the CCD). 

Numerically, cr v from this calculation is 16 m/s (with d = .58, Iq = 10 5 ). The finite 
spectral resolution, however, tends to wash out the signal. For Gaussian transmission- 
profiles (fwhm = 100 mA), a v is 25 m/s. For the transmission-profiles in figure 2.5, 
cr v is 30 m/s as shown in figure 2.14 (a-only method) and figure 2.15 (a, 3 method). 
The scatter in figure 2.15 is due to the method by which the calibration table 3 ) 


is calculated. 

The sudden increase of cr v for velocities above 3500 m/s is due to the increasing 
dependence of t* on 3 together with the fact that 3 is only allotted 9 bits in the 
calibration-table address. Allotting more bits to 3 would proportionately reduce this 
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Figure 2.15: Statistical error a v vs. velocity for different line-widths (a, /? method). 





CHAPTER 2. MEASUREMENT TECHNIQUES 


29 


behavior. 


2.2.5 D iscussion of the Method 

The studies indicate that velocity errors due to line-depth variations and Lyot filter 
drift can be minimized by proper design of the Lyot filter profile. In principle it 
should be possible to eliminate the errors entirely, though only a crude optimization 
has been done here. 

Line-width variations can produce velocity errors of up to a few hundred m/s 
when the velocity is calculated trom a single rational function of the data (or-only 
method). The error is reduced to a few m/s by using two rational functions (a, /? 
method). The more accurate computation requires an extra image buffer and entails 
doing 2.2 times as many operations. 

A line-asymmetry will affect the measured velocity because the line-bisector is 
sampled at different depths across the solar disk. This happens because the average 
velocity is a function of disk position, and the data tend to sample more of the 
line-wings for large Doppler shifts. 

Data taken at La Palma in 11)88 show that the averaging of subresolution velocity 
elements produces a ’‘noise'' level of - ±20 m/s. Hence, it seems reasonable to set 
this as an upper limit for systematic errors. 

The statistical error for a four point measurement is ~ ±30 m/s with a CCD 
well-depth of 2 x 10 5 photons. The error becomes significantly larger at v ~ 3500 
m/s for the a, 3 method. This can be rectified by use of a larger calibration table. 
The errors for 2 pass (8 point) and 4 pass (16 point) measurements are 21 m/s and 
15 m/s, respectively. 

In order to minimize t he statistical error, sidelobes should be kept to a minimum. 
Hence the optimum Lyot filter-profile would have a flat peak with wings just large 
enough to ensure integration of equal amounts of continuum at each wavelength 
setting. 

This recommendation tor the Lyot filter width has been used in the design of the 
MDI instrument. The choice of velocity calculation method is still pending further 
study. It is not yet* known whether line-width or line-asymmetry variations will 
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contribute more to measurement errors. Only it line- asymmetry variations are small 
can the a, 3 method used. 


2.3 Magnetograms 


2.3.1 Inverse Zeeman Effect 


Spectral lines of atoms in a magnetic field are split by the Zeeman effect. For a field 
component H\\ pointing toward the observer, an emission line of wavelength A 0 is split 
into two components (a components) with opposite circular polarizations. The left 
and right circular polarizations (LCP arid RCP) are displaced from A 0 by — A,\ and 
T A \ , respectively. 


A a 




4 irnifC 2 


A 0 2 (g'M'-g"M“ 


(2.17) 


where g' and g" are the Lande factors, and M' and M" are the magnetic quantum 
numbers for the initial and final states (Bray k Loughhead 1964). The magnetic 
quantum numbers are subject to the selection rule M" — M' — ±1. For the Fe I 
(6302.508 A) line. 


A\ 

"ll 


47.4 mA/k 


gauss. 


(2.18) 


For absorption lines, the LCP and RCP components will be preferentially absorbed 
at A 0 — A,\ and Aq + A,\. respectively. This is termed the inverse Zeeman effect. 


2.3.2 Longitudinal Field Measurement 


From figure 2.16 we can write down the measured LCP and RCP intensities. 

Id/ 1 


h = lo + 

//?=/( 0 — 


d\ 

di 


a a 


dX 


(2.19) 

( 2 . 20 ) 


where is assumed to be constant over the range \ w ± A,\. 
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the Fourier transform. Consider a monochromatic plane wave p with wave vector 
k 0 = ( ko x , k 0y ) and frequency ^o: 

f){xj) - cos {h'Q • X - uJ 0 1) (2.24) 

The Fourier transform of p is defined as, 

p } = J f J cos (4 ‘ x - ^ 0 t)e~ t ^ x+UJt) d 2 x dt (2.25) 

Evaluating the integral, 

{/*} = ^ — A‘q ) <^(u; — u^o) T + k'o) 6(uj + u,’o) (2.26) 

The two delta functions, lying opposite one-another in (fc, a;) space, illustrate a general 
property of Fourier transforms. For real valued functions, the (complex) Fourier 
transform will have a factor of 2 redundancy in information. In fact, we only need 
to look at a half-space (containing the origin) to obtain all the information about 
the original function. It is conceptually easiest to eliminate the negative frequency 
half- space, and consider only: 

— DC < h\ r < DC, — DC < k y < DC, 0 < u; < DO 

It is important to account for relative motion of the observers field of view with 
respect to the observed medium. If the medium is moving at constant velocity v with 
respect to the observer, we have x — > x — vt. This produces a measured frequency 
shift jj — ► u? + for the plane wave, where 

= (2.27) 

A final and extensively used property of Fourier transforms is the convolution 
theorem. It states that the Fourier transform of the convolution of two functions 
f(x) and g(x) is equal to the product of the Fourier transforms of the functions, 

/(0</U - 0^} = Hf(x)} F{g(x)} (2.28) 
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2.4.1 Practical Considerations 

The discrete Fourier transform, or Fourier series, can be viewed as an extreme example 
of a Fourier transform (Bracewell 1986). The series representation Fk is useful for 
discrete, evenly spaced data f n with periodic boundaries. It is defined by, 

Fk = F {.fn} = f n e~ 2i * nk ' N (2.29) 

V n = l 

The inverse transform is defined by, 

.fn = F~ l {F,} = ^'t F k e 2 '™ k ' N (2.30) 

* v k= i 

For measurements made at an interval <*>, the maximum Fourier frequency that can 
be represented is the Nyquist frequency. 



The discrete transform has the advantage that it can be evaluated using a very 
efficient algorithm, the Fast Fourier Transform (FFT) algorithm. The data from the 
series of Doppler and intensity images is indeed binned discretely on an evenly spaced 
grid in the spatial dimensions (by pixel) and temporal dimension (by image number). 
The condition of periodic boundaries must be enforced, however, by apodizing (see 
below). 

There are several algorithms available for discrete Fourier and related transforms. 
For handling real valued data, it is more efficient (by a factor of 2) in computing 
time and memory usage to construct the Fourier transform from a Fast Sine-Cosine 
Transform. For doing convolutions to smooth and filter data, it is most efficient and 
conceptually simplest to use the Fast Hartley Transform. The abbreviation FFT is 
used generically hereafter to refer to the Sine-Cosine and Hartley transforms. 

The FFT algorithms used were designed for data series containing N points, where 
N is even and factorable into low primes. The spatial dimensions were 512 pixels in 
both x and y ( N x = N y = 2 <J ). The temporal dimension was 480 images {N t = 2 5 • 3 ■ 5 ) . 
This number was chosen because it entailed the largest “nice” number of filtergrams 
that would fit on a single data tape. 
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2.4.2 Apodization 


To approximate periodic boundaries, the data is multiplied in each dimension by a 
function which tapers to zero at. the edges. The exact form of the function must 
be chosen very carefully to avoid disturbing the Fourier spectrum. The convolution 
theorem tells us that multiplying the data by an apodizing function is equivalent to 
convolving the Fourier spectrum with the Fourier transform of the apodizing function. 
Therefore, the apodizing function should be constructed such that its transform has 
a narrow bandwidth and very small sidelobes. 

Several standard apodizing functions are listed in Harris (1978). The 3-term 
Blackman- Harris function was chosen. Its discrete Fourier transform has a 3 dB 
width of 1.56 bins, and largest sidelobes of —61 dB. The windowing function for N 
bins is given by 


IF,, = flu + <i\ cos 


y)+« 2 cos(f„) 


(2.32) 


lor n — — ^ T 1 * . • • i — 1 

where r/ 0 ~ 0.44959 
ai = 0.49364 
a 2 = 0.05677 

When filtering data as done in chapter 7, it is not, in principle, necessary to 
apodize the data. For the La Palma data sets, however, the image rotates during 
the observations. Removal of t his rotation from the images leaves the corners of the 
images without real information. To minimize round off errors from the corners, the 
data were masked before filtering. The mask used had a cosine taper (to zero) at the 
edges of a flat central disk: 


— 


cos 


1 if r mn < N - 32 

0 if r mn > N 

( — (r mti - N + 32)) otherwise 


(2.33) 


where r 


= \/m 2 + 77 2 , 


A V i 

= —J T — 1 
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2.4.3 Normalization 


The FFT’s are normalized such that the spectral power of a sine wave is equal to the 
square of the wave-height before apodization. 


T 2 

Y t sin A’o” } | =1 

k=-$ 


(2.34) 


The exact normalization is a weak function ot k and is different for the real and 
imaginary parts of the transform. However, it can be approximated by a constant. 


= -MF{Wnfn}\ 2 (2.35) 

(2.36) 

This approximation is valid tor wavelengths that are short compared to the scale 
of variation ot the apodizing (unction H as is the case for all wavelengths and 
frequencies of interest. 

There is an additional wavenumber-dependent normalization factor that corrects 
for the attenuation ot spatial trequencies by the atmosphere and optics. This correc- 
tion is described in sections 3.4 and 5.4.1. 
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2.4.4 Statistical Error Estimates 


An error estimate can be obtained lor the spectral power |E.| 2 by a standard propa- 
gation of error calculation. II we assume a constant error oy on the data / n , we can 
use equation 2.29 to find the error on the power a |F fc p- 


a \F k V 



(2.37) 


In three dimensions, this becomes: 


°y 


If*. I 




(2.38) 



I i 



Chapter 3 

Telescope and Instrumentation 


3.1 Telescope 

The 50 cm Vacuum Solar Telescope at the Swedish Solar Observatory is shown 
schematically in figure 3.1. A photograph ot the tower at Roque De Los Mucha- 
chos on the island of San Miguel De La Palma (Canary Islands, Spain) appears in 
figure 3.2. The telescope is a refractor with an / number of 45. 

The image forming lens is a 50 cm achromatic doublet. Two flat mirrors Ml and 
M2 in the altazimuth turret guide the beam down an evacuated tube. The first lens 
of the achromat lorrns the vacuum window. At the bottom, the flat mirror M3 directs 
the beam through a vacuum window to the (Ursies) observing room. Extreme care 
was taken by the designers of the optics to ensure the telescope would perform at the 
diffraction limit (Scharmer et nl. 1085). 

3.2 Postfocus Instrumentation 

3.2.1 Tunable Filter 

Lyot Filter Operation 

A Lyot filter is constructed from a series of simple birefringent elements. A represen- 
tation of one such element is shown in figure 3.3. An incident pulse of light is polarized 




Figure 3.1: Schematic diagram of the Swedish 50 cm Vacuum Solar Telescope. From 
Scharmer et ai (19S5). 

at 45° to the optical axes of a birefringent crystal. Fast and slow pulses emerge from 
the crystal, and pass through an exit polarizer. The output of the element is, thus, 
a pair of identically polarized, equal amplitude pulses, with a temporal separation 
proportional to the crystal thickness. The spectral response of such an element is 
given by the Fourier transform of its impulse response. The Fourier transform of a 
pair of delta functions is a cosine. 

A Lyot filter consists of a cascade of elements, each double the length of the 
previous one (see figure 3.4). The transmission of each element as a function of 
frequency is shown in figure 3.5. The longest element has the most closely spaced 
maxima. The next longest has minima that fall on every second maximum of the 
longest element, and so forth. The net effect is to produce a transmission profile 
with a mainlobe width approximately equal to that for the longest element; and a 
free spectral range, or distance between peaks, equal to the period for the shortest 
element. 
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Figure 3.2: Swedish Solar Observatory. La Palma, Spain. 
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Figure 3.3: Pulse propagation through a simple element. From Title & Rosenberg 
(1981)1 1 



Figure 3.4: Lyot filter with four elements. From Title & Rosenberg (1981). 
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Figure 3.5: Transmission vs. trequency for the components of a Lyot filter (A, B, 
C, D) and the transmission through successive stages of the cascade (E, F, G). From 
Title & Rosenberg (1981). 
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Figure 3.6: Birefringent tuning element. From Title h Rosenberg (1981). 


Tuning a Birefringent Element 

The bandpass of a birefringent filter can be tuned over the full free spectral range 
by individually tuning the component elements. The best practical method for doing 
this is illustrated in figure 3.6 (Title Rosenberg 1981). As in figure 3.3, a fast 
and slow pulse emerge from the crystal parallel to the fast and slow axes. The delay 
between the pulses is inversely proportional to wavelength. Reconstituting the two 
polarizations gives elliptically polarized light. The axes of the ellipse are horizontal 
and vertical, and the ratio of the axes, the ellipticity, is a function of wavelength. 

A quarter wave plate transforms the polarization from elliptical to linear, and the 
orientation of the resulting polarization vector is dependent on ellipticity. Finally, 
a rotatable half wave plate and fixed exit polarizer are used to select the linear 
polarization, and hence, the wavelength. 
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The Solar Optical Universal Polarimeter 

The Solar Optical Universal Polarimeter (SOUP), developed by LPARL, employs a 
nine element tunable birefringent filter. It is similar in design to a Lyot filter, but each 
birefringent element contains a partial polarizer to apoclize its transfer function. The 
result is a filter with smaller sidelobes than a Lyot filter. The SOUP instrument was 
used on the Spacelab 2 shuttle flight, and the filter has become the evaluation model 
for the Orbiting Solar Laboratory — Coordinated Instrument Package (OSL-CIP). 
Software-selectable blocking filters, each 6-8 A wide, are centered on solar absorption 
lines. Waveplates can also be selected to measure left and right circular polarizations 
as well as linear polarizations spaced every 45°. 

At a wavelength of 5896 A. the free spectral range of the tunable filter is 14 A, 
and the f whin of the passband is 62 mA. A "wide mode is also available by rotating a 
half wave plate within the longest element. This effectively shortens the element by a 
factor of 3, and increases t he bandwidth to 103 mA. The stepping motors that rotate 
the tuning half wave plates permit tuning in 3.14 mA steps. Tuning is repeatable to 
within half a step. 

Temperature gradients and imperfections in the filter components produce small, 
but measurable, variations in filter characteristics over the field of view. The center 
of the passband varies by up to 10 mA and the Jivhm, by up to 5 mA over the field of 
view. In addition, the total amount of light transmitted (integrated over wavelength) 
is modulated as the filter is tuned. (See figure 3.8.) This modulation accounts for 
most of the error in the velocity and intensity measurements. 


3.2.2 CCD Camera 

The CCD (charge coupled device) camera uses a Texas Instruments 1024 x 1024 pixel, 
virtual phase detector. It is also an evaluation model for OSL-CIP. Subareas of the 
full detector can be read out to reduce the time required to write data to tape. In 
addition, 2 x 2 groups of pixels can be summed to decrease exposure time by a factor 
of 4. This 2 x 2 summing mode was used for the majority of observations listed in 
chapter 4. 
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Figure 3.7: Transmission profile of the SOUP tunable filter vs. wavelength. Trans- 
mission was measured using a spectrometer, with the filter tuned to 6328 A. 

The detector geometry consists of a rectangular grid of pixels. Each pixel is 
efficient at collecting light up to the boundary with its adjacent pixels. Thus, the full 
area of the CCD array is sensitive to light. 

The charge accumulated on a pixel is proportional to the integrated light falling 
on that pixel. The signal also integrates a current present in the absence of light, 
the w dark current.” The dark current depends on device temperature. To keep this 
effect constant, and to a minimum, the CCD is kept in a vessel at low pressure, and 
is cooled thermoelectricallv to —20° C. 

3,2.3 Image Stabilization System 

Compensation is made for small image displacements, caused by seeing and inaccurate 
telescope tracking, by an active tilt mirror in the beam path. The mirror is moved 
in accordance with signals from a sunspot tracker which senses image motion. The 
active mirror is mounted on three piezoelectric actuators which can tilt the mirror on 
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Figure 3.8: Relative transmission of the SOUP tunable filter vs. tuned wavelength. 
Transmission was measured using a continuum illumination source. 

two axes. The system has a bandwidth of 65 Hz. 

The spot tracker is an array of four photodiodes (a quad-cell). A sunspot image 
centered on the quad-cell will produce equal signals from each diode. If the image 
moves toward one of the four quadrants, however, that quadrant will experience less 
illumination and produce a smaller signal. The differences of signals from the four 
quadrants are used to determine compensatory signals for the active mirror. The 
integrated compensatory signals are used as a feedback to the telescope tracking 
system. 


3.3 Optical Layout 

The optical bench was set up in the Ursies room as shown in the photograph in 
figure 3.9. The schematic layout of components is shown in figure 3.10. The beam 
from the telescope strikes the active mirror, and is deflected through a field stop 
slightly larger than the active area of the CCD. The heat stop, used to filter out 
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infrared wavelengths, was removed from the light path because it produced fringing 
patterns. 

Wavelengths of 5000 7000 A are reflected by the dichroic beam splitter to the 
tunable filter path. Infrared wavelengths are not reflected. Other wavelengths are 
reimaged through the (’a II K-line filter and the 4308 A bandpass filter to the K-line 
and “white light" cameras, respectively. The white light images were used to monitor 
the seeing conditions and the focus in real-time. The K-line camera was not used for 
these observations. 

After reflection from the dichroic beam splitter, the beam is collimated by a lens 
on a motor driven stage. The magnification of the system is chosen bv substituting 
lenses of different focal length. The stage is moved automatically to compensate for 
the chromatic dispersion of the telescope, and the whole assembly, moved manually to 
adjust the focus. The beam passes through an aperture stop and is reimaged through 
a set of selectable waveplates and blocking filters, and the tunable filter. For solar 
observations, the beam passes through a sector shutter to the CCD camera. To check 
the position of a spectral line, t he flip mirror was inserted into the beam to deflect it 
to the photometer. 

Part of the beam is diverted, before the polarizer wheel, to the spot tracker. There 
is a neutral density filter and adjustable zoom lens system in this path to optimize 
the intensity and magnification of the image for the spot tracker. 

The motor driven stage, polarizer wheel, blocking filter wheel, and tunable filter 
were controlled by a (PC compatible) COMPAQ computer. The shutter and CCD 
camera were controlled by a VAX Station 3200. The system was free-running, with no 
timebase to control the data acquisition. The timing of the images, therefore, varied 
as shown in figure 3.11. The time mark given to each image came from a crystal 
controlled clock, and is accurate to — 1 ms. 
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Figure :}.9: Optical bench setup. 
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Figure 3.11: Timing offset for images vs. time from start of run. The offset is the 
difference between the image time mark and a straight line fit to the time marks. 

3.4 Modulation Transfer Function of the System 



The modulation transfer function (MTF) of an optical system describes the effect of 
aberrations and diffraction on an image formed through the system. If I g {x) is the 
image intensity distribution expected from aberration-free geometrical optics, then 
the real image intensity distribution / r (.c) is given by, 

M. f) = .F -1 {MTF(£) • F{I g {x)}} (3.1) 

The attenuation of the spectral power |j r {/ 3 }| 2 is given by the square of the MTF 
functions. 

The components which contributed significantly to the total MTF are the entrance 
pupil, the detector pixel shape, and the earth’s atmosphere. The total MTF is the 
product of the MTF’s of these individual components. The atmospheric MTF will be 
treated in section 5.4.1. 
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3.4.1 Entrance Pupil 

The entrance pupil is defined by the diameter d (48 cm) of the heat shield diaphragm 
in front of the objective. Idle MTF of the pupil is given by the autocorrelation of the 
pupil function (a disk). For wavenumbers, 0 < k < A' Nvq , 

MTF(A-) = - [cos -1 k - Ky/l - k 1 2 ] (3.2) 

A / # pixels \ k 

where k = — - — — — ; — 

2d \ field oi viewy k^ yq 


3*4.2 CCD Camera 

The effect of summing light over the area of a CCD pixel is equivalent to convolving 
the image distribution with the pixel shape, and sampling the new distribution at 
points centered on each pixel. Thus, the MTF of the detector is the modulus of the 
Fourier transform of the pixel function (a square). 


MTF( k) = sine 


sine K y 


(3.3) 


where k = 


1 k 

^ ^Nyq 


SU17TX 

and sinc.r = 

7TX 

This can be approximated by a cylindrically symmetric function with no more than 
a 2% error. 


MTF (k) — sine k 


(3.4) 


The square of the MTF functions are plotted in figure 3.12 for the entrance pupil and 
CCD detector. 


MTF 
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Figure 3.12: Squared MTF of telescope and CCD detector 
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Chapter 4 
Observations 


The original topic ol this dissertation was “sunspot seismology, that is, the study 
waves propagating into and out ot sunspots. That topic was set aside when waves were 
discovered at substantially higher wavenumber and frequency than had previously 
been observed. The new topic of analysis became the characterization of properties 
of these waves. Man)' ot t he observing runs are, therefore, extraneous to the analysis 
described here. They have been included in the list of data in this chapter, but no 
calibration or reduction has been performed on them. 


4.1 Observing Procedure 

The observing runs consist ot a repeating sequence of filtergrams with a temporal 
spacing of about 3.5 s between filtergrams. The wavelength, filter width, polariza- 
tion, and relative exposure time for each filtergram in the sequence is defined in the 
sequence file residing on the COMPAQ computer that runs the filter. 

Calibration data was taken before and after observing runs. The order of calibra- 
tion and observations was as follows: 

1. Dark Currents: 4 measurements of dark currents for each exposure setting. 
Done at the beginning and end ot each day. 
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2. Line Centering: Definition of the nominal line center from which wave- 
length offsets are measured. The photometer detector was used to inte- 
grate light over the field of view with the target region centered in the 
field. 

3. Slide Projector Flats: 4 repetitions of flat field measurements using slide 
projector source. 

4. Observing Run: One or more tapes-worth of observations of target region. 

5. Solar Flats: 4 repetitions of flat field measurements using defocused sun- 
light source. 

6. Slide Projector l lats: Repeat of step 3. 

The line centering step is performed by scanning the filter through a range of 
wavelengths (from short lo long), and plotting intensity versus wavelength. This step 
was not performed correctly due to a design flaw in the photometer interface with 
the COMPAQ computer. The computer did not wait for the full signal integration 
time before reading the light intensity. In the resulting line profile scan, the intensity 
points lagged behind the wavelength points. This caused the line to appear shifted 
to longer wavelength. The line' was, therefore, “centered 11 to an erroneously short 
wavelength. This offset of — 10 mA for the Na D\ line was discovered by making a 
line scan using the CCD camera instead of the photometer (see figure 4.1). 


4.2 Calibration Method 

The purpose of calibration is to produce images with a signal to illumination ratio 
that is independent of both the pixel location and the wavelength setting of the 
tunable filter. The procedure described here is more complicated than the standard 
one used by the Lockheed group. This is necessary to provide an accurate calibration 
for dopplergrams. The flat lidding procedure calibrates the whole system, optics and 
detector together. 

The CCD is a very linear device, provided the light intensity stays below saturation 
levels. Response of a pixd can be characterized by a gain and an offset. The offset 
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Figure 4.1: Line scan of Xa Di Line using the CCD camera. Signal, integrated over 
the field of view, vs. wavelength offset. 

is determined by the amount of charge that accumulates on the CCD in the absence 
of illumination (the Mark current" ). It must be measured for each exposure time. 
Each pixel has a characteristic gain and dark current. The tunable filter is also 
responsible for some variation of the signal strength as a function of wavelength A 
and spatial position x. The net response of the system is to produce filtergrams F % (x) 
at wavelengths A t1 given a source intensity distribution S(x, A). 

F,(x) = jS{x,\)A,{x,\)d\ + D{x) (4.1) 

where Ai(x, A) is the product of filter transmission and CCD gain at filter setting A,, 
and D(x) is the CCD dark current. 

Because the filter bandpass is narrow compared with spectral features in the 
source, we can approximate: 

J S(.v, A) Ai(x, A) dX = 5(x, A,) J Ai(x,\)d\. (4.2) 



Thus 


Fj(x) = S{x, A,) BJx) + D(x) 


(4.3) 
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where B,(x) = jA t {x,\)d\. (4.4) 

In order to determine 5\ the quantity of interest, we must first measure B x and D. 
The dark current D is easily measured by making filtergrams with no illumination 
source. For 5 = 0, we have D = F t . 

If we had an illumination source that were both flat in spectrum and flat across 
the field of view, we could measure B t directly. The sources available, however, were: 

1. Defocused Sunlight 5',sdA): Flat across field, but not spectrally flat. 

2. Slide Projector Lamp S’/fi.r): Flat in spectrum, but not flat across field. 


The technique is to make one filtergram Ff at wavelength A c in the continuum 
(away from the absorption line), using defocused sunlight. And, to make filtergrams 
F C L at A c and F t L at A; usi ng the slide projector lamp. This gives us the following 
three measurements: 


E s (.v) 

= S s (\ c )B c (.7) + D s (x) 

(4-5) 

E L (T) 

= S L {?) B C (T) + D l (x) 

(4.6) 


— ^ l ( x ) B t (x) + Di(x) 

(4.7) 


Given the dark currents /4/, and Ds, we can solve for the gain correction Q t at A, 
defined as Cy, ( .r) = Bi(x)~ l . 


C7,(.P 


[ft ~ Dl] 

[ F l _ Dl][f s _ Ds] 


(4.8) 


There is some arbitrary overall normalization which depends on 5s(A c ). In practice 
the spatial average of Q x at the \ c filter position was defined as unity, and the other 
gain corrections, normalized accordingly. 

With the gain corrections and dark currents in hand, we can determine arbitrary 
source distributions from filtergrams. 


>'(.r. A,) = C, t (.r) [F,(.r) - D(x)\ 


(4.9) 
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4.3 Solar Observations 

The sequence of filtergrams for the sunspot seismology study was named SEISNAFE. 
It was designed to determine the following quantities: 

1. Doppler shift in mid-chromosphere (in sunspot umbrae). 

[H I (Balmer) n — line core] 

2. Doppler shift in photosphere. 

[Na. D ( — wings] 

3. Doppler shift in low chromosphere. 

[Na D | — core] 

-f. Doppler shift in photosphere. 

[Fe 1 (5576 A)] 

5. Continuum intensity in green and red. 

[5575 A, 589 1 A] 

6. Photospheric magnetic field (line of sight component). 

[Fe 1 (6 302 A)] 

The sequence repeated every 17 s. Filter parameters appear in table 4.1. T he corre- 
sponding calibration sequences were NAI EFLAT for solar flats, and SEISSP for slide 
projector flats. The observing runs of appreciable length are listed in table 4.2. 

The sequence of filtergrams used for the high-frequency oscillation study was 
named NAFAST. It was designed to make dopplergrams in the low chromosphere 
and photophere, using filtergrams in the core and wings of Na Di, respectively. All 
other measurements were eliminated to increase the repetition rate to once every 
17.7 s. The sequence parameters appear in table 4.3 and the corresponding calibra- 
tion sequences were NAFF FLAT for solar flats, and NAFASTSP for slide projector 
flats. 

The NAFAST sequence was run only once on August 2, 1990 beginning 13:51 
UT. The target region was an active region (AR61S1) containing a small, decaying 
sunspot. This is shown on a full-disk magnetogram in figure 4.2. The run lasted 
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Table 4.1: Filter ' 

Parameters for SEISNAFE Sequence. 

Line (A) 

Offset 

(niA) 

Filter Mode 
(width) 

Polarization 
(0 = linear) 

Purpose 

H I a 6562.8 

-100 

0 

-100 

wide 

0 

dopplergram 

Na Di 5895.9 

-to 

0 

10 

narrow 

0 

doppiergram 

Na Dj 5895.9 


narrow 

0 

dopplergram 

Na Di 5S95.9 

-2000 

narrow 

0 

continuum 

Fe I 5576.1 


wide 

0 

dopplergram 

Fe I 5576.1 

-750 

wide 

0 

continuum 

Fe I 6302.5 


narrow 

LCP 

RCP 

magnetogram 














CHAPTER 4. OBSERVATIONS 57 


Table 4.2: La Palma Observing Runs using SEISNAFE Sequence. 


Date 

(1990) 

Time 

(UT) 

Du ration 
(min) 

Seeing 

FOV 

Target Region 

28 Jul 


■B 

good/v.good 

75" 

AR6168 (medium sunspot) 



m 

fair/ var. 

150" 

AR6168 

29 Jul 

8:39 

156 

fair/var. 

150" 

AR6174 (medium sunspot) 


16:53 

126 

good/ var. 

150" 


30 Jul 

7:50 

281 

variable 

224" 



13:12 

293 

variable 

224" 


31 Jul 

8:34 

237 

good /var. 

224" 

— flare in morning 


13:09 

347 

variable 

224" 


1 Aug 

8:08 

256 

good /var. 

224" 

AR6181 (small sunspot) 

1 Aug 

Egga 

310 



AR6179 (big sunspot) 

2 Aug 

mm 

151 



— overexpd for umbral osc. 

\mwxm 

10:52 

101 

fair/ var. 

224" 

AR61S6 (big sunspot on limb) 

5 Aug 

9:20 

190 



quiet region near AR6192 
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Figure 4.2: Full-disk magnet ogram from the National Solar Observatory at Kitt Peak, 
The field of view containing ARG181 is outlined. 

143 min with good/ variable seeing that deteriorated in the last 30 min of the run. To 
obtain more information about t lit' active region, the SEISNAFE sequence was run 
before and after the XAFAS 1 runs. A continuum image is shown in figure 4.3, a high 
resolution magnetogram, in figure 1.1. and, an offhand H o image, in figure 4.5. 
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Figure 4.3: Continuum intensity (5575A) image of AR61S1. The small tick marks 



The small tick marks are 5". 
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Figure 4.5: H a image of AR61S1 (400 mA in the blue wing). The small tick marks 



Chapter 5 

Data Reduction 


The raw data used for this dissertation reside on a (full) 2.2 Gigabyte Exabyte tape. In 
producing a power spectrum, the volume of data is reduced by a factor of ~ 3 x 10 4 . 
There is more to the reduction procedure, however, than extracting a manageable 
amount of information from a cumbersome data set. The quality of the data depends 
on careful calibration to account for instrument response, image rotation and drift 
caused by telescope motion, and, distortion caused by the atmosphere. The absolute 
spatial scale of the images must also be determined accurately. This chapter describes 
the details of and cross-checks for the data reduction procedures. 

The reduction steps are listed below in the order they were performed: 

1. Flat Fielding: Calibrate for variation in filter transmission and CCD gain, 
and subtract CCD dark current. 

2. CCD Flaw Correction: Interpolate over nonfunctional pixels. 

3. Correction for Atmospheric Blurring: Normalization of measured MTF. 

4. Rigid Rotation: Compensate for image rotation caused by telescope track- 
ing motions during observations. 

5. Correction for Atmospheric CShimmerC (Destretching) 

6. Calculation of Dopplergrams and Line Center Intensity. 

7. Rigid Alignment: Compensate (or image drift caused by telescope tracking 
errors. 
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S. Fourier Transform of Dopplergrams and Intensities: Described in 
section 2.4. 

9. Angular Scale Determination: Convert units from pixels to arcseconds. 

10. Distance Scale Determination: Measurement of foreshortening and direc- 
tion toward disk center. Described in section 6 . 1 . 1 . 

11 . Measurement of Oscillation Properties: Described in chapters 6 and 7 . 

5.1 Flat Fielding 

The flat fielding scheme is explained in section 4.2. We can use equation 4.9 to 
calibrate filtergrams, given dark currents D and gain corrections Q{. The dark current 
is an average of lour filtergrams made with a. light stop placed in front of the filter. 
The gain correction is constructed using equation 4.8 as follows. 

The three flat field measurements Ff , Ff , and F t L are each made from an average 
of four filtergrams. Because the filter stepping motors did not always find the correct 
positions, it was sometimes necessary to throw out one of the four measurements 
and use only three to make the average. A filtergram was thrown out if the spatial 
average of the central 64 x 64 pixels of the image differed by more than 5 cr from the 
spatial averages of subareas of the other three images, where a is defined by inefficient 
statistics. 

■ 7 , = W ~ /min (5.1) 

N i 

The gain correction is the product of (F c 5 — F 5)" 1 and (F C L — Di)j{F ) L — Di ). 
The latter is the ratio of the transmittance in continuum to transmittance an X t . 
This ratio should be slowly varying over the field of view. The dark currents are 
not constant with time, however, and the very long exposure times required for slide 
projector flats (2-3 s) caused accumulation of large charges at some small groups of 
pixels. It was, thus, necessary to “deflaw" F^ and F t L (see next section), and smooth 
the transmittance ratio using a Gaussian point-spread function of 18 pixel fwhm. 

The procedure described here produced (flat fielded) images with signal to illu- 
mination ratios that were constant to within about 3 % across the field of view, and, 
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more importantly, to within about 5% across the wavelength range of an absorption 
line. These numbers were measured from filtergrams of defocused sunlight (the flat 
field illumination source) at continuum wavelengths. 


5.2 CCD Flaw Correction 

Of the 2.5 x 10 5 pixels (in 2 x 2 summing mode), 4513 were not functional. An 
additional 491 pixels had exceptionally high dark currents for exposure times greater 
than about 500 ms. To minimize problems caused by the missing data, values were 
filled in at these pixel locations by linear interpolation from adjacent pixels. This 
procedure is termed "deflawing . 

A good map of nonfunctional pixel locations was available for exposure times 
below 500 ms. For the longer exposure times, a map was constructed using the dark 
current runs. To do this, a cubic trend was fitted to each row of a filtergram. Pixels 
which differed from the trend by more than 7% were flagged as “bad”, and their 
locations, entered in the map. 

5.3 Image Rotation 

The altitude-azimuth axes of mirrors Ml and M2 cause the solar image to rotate 
during the day as the telescope tracks the sun. The rotational direction and rate are 
functions of date and time of day, and can be calculated from geometrical consider- 
ations. Standard software exists to determine the angle of rotation and regrid the 
images by bicubic interpolation. This removes the apparent rotation in a sequence of 
images, but some data is lost from the corners of the images. 


5.4 Atmospheric Distortion 

Turbulence in the earth's atmosphere will modify wavefronts of light before they reach 
the telescope. Any plane wave passing through the atmosphere will be distorted. 
The resulting wavefront, however, can be approximated by a set of planar patches. 
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The isoplanatic angle describes the angular extent over which image distortions are 
correlated. For excellent seeing conditions, it can be over 100" (von der Liihe 1988). 
For poor seeing, it can be less than 1". 

A change in angle of incidence of a planar wavefront produces an image displace- 
ment. When different, parts of an image are displaced in different directions, the 
image appears distorted. The distortions change with the atmospheric turbulence. 
This can occur over a range of timescales, from 1 to 1000 ms, with most power in the 
10-100 ms range. 

5.4.1 MTF Normalization 

The distortion of a wavefront from a given object point differs over the telescope 
aperture. Moreover, local displacements of up to an arcsecond cause the image to 
move around in times of order 10 ms. For exposures lasting tens to hundreds of 
milliseconds, these changing distortions will be superimposed to produce a blurred 
filtergram. The effect is equivalent to convolving the image with some unknown 
point spread function. The Fourier transform of the point spread function is called 
the MTF. It is assumed here that the point spread function is the same everywhere 
in the field of view of a given filtergram. 

The atmospheric MTF can be estimated from the spatial frequency spectrum of 
a filtergram S(f). A contour plot of spatial frequency spectra is shown in figure 5.1 
for the —40 mA filtergram set. The sharpest filtergrams in the set are assumed to 
be diffraction limited (atmospheric MTF = 1). Since the structure of solar images 
should not change appreciably over the observing run, the spatial spectra of diffraction 
limited images should be a constant S max {f). Any difference between the spectrum 
of a given filtergram and 5 max must, therefore, be due to the atmospheric MTF. 

S(f) = MTF alill (/) ■ 5, nax (/) (5.2) 

In practice the log-spectra were compared with the mean log-spectrum log S. 
Log-spectra of individual filtergrams were smoothed by convolution with a cosine-bell, 
and the difference between log S and t he smoothed spectrum was used as a correction 


Log(RMS amplitude) 
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Figure 5.1: Spatial frequency log-spectra of Na Dj -40 mA filtergrams. The solid 
central line is the mean log-power; the dashed lines are lcr contours; the outer dotted 
lines are the maximum and minimum power at a given frequency of any image. 
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factor to the filtergram frequencies. Thus, the ”MTF normalization” procedure forced 
the spatial spectrum of each filtergram to look like the average spectrum. 

In principle the individual filtergrams could have been normalized to 5 max . This 
would have involved very large correction factors for the poorer filtergrams, and would 
have boosted the noise in t hese filtergrams considerably. Noise in filtergrams produces 
noise in dopplergrams and line center intensities in a nonlinear manner. So, instead, 
the difference between the mean spectrum and 5 max was retained as a correction 
factor to be applied to the oscillation power. (See chapter 6.) As a final caveat, it 
should be noted that even the best images, assumed above to be diffraction limited, 
actually include some atmospheric blurring. Hence, the correction factor will provide 
an underestimate for power at the highest wavenumbers. 

The MTF normalization procedure was checked by comparing filtergrams taken 
in poor seeing with filtergrams taken at almost the same time in good seeing. After 
normalization, both filtergrams showed the same solar structures. 

5.4.2 Dest retching 

The displacements caused by atmospheric turbulence vary over the field of view. 
They are independent of one another over angular scales greater than the isoplanatic 
angle. Relative displacements between successive images can be measured by local 
correlation tracking (November 19S6). This procedure finds the displacement of a 
small subarea of one image with respect to the other from the cross-correlation of the 
two. A displacement vector field is measured using a patchwork of subareas of the 
images. A relative displacement grid of 64 x 64 vectors was made, correlating each 
filtergram with the preceding one. This was done using standard software. 

Local displacements are present not only because of seeing, but also because of 

wave motion and horizontal flows on the sun. We wish to remove the seeing displace- 
ments, but keep the real solar motions. It is possible to distinguish the two on the 

basis of timescale of the motion. Horizontal flows have characteristic timescales of 

hours, oscillations have timescales of minutes, and seeing has a timescale of tens of 
milliseconds. For filtergrams taken seconds apart, seeing motions will appear random, 
while solar motions will be coherent. 
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The seeing displacement grid was determined from the relative displacement grid 
as follows. For each vector component at each grid location, the relative displacements 
of preceding filtergrams were summed to produce a total net displacement (from the 
first frame). The total displacement was smoothed by convolution with a Gaussian 
of 90 s (5 bin) fwhm , and this smooth trend, subtracted to produce the seeing grid. 

It is important to have high temporal resolution in the data. Correlation tracking 
can lock onto the crests of waves. So, if the impulse response function, used to smooth 
the data, is wider than ~ 120 s, the oscillations will be removed from the data along 
with the seeing. The standard Lockheed method for finding the seeing grid is designed 
for horizontal flow studies. It uses boxcar averages of 20 or more bins. This would 
be inappropriate for seismology studies. 

With the seeing displacements measured, the images can be “destretched” to 
their undistorted states. The algorithm to do this treats images as rubber sheets, and 
stretches the sheets (by regridding) based on the local displacement vectors. This part 
of the process also used standard software. The set of filtergrams for each wavelength 
was destretched independently of the other wavelength sets. Destretched filtergrams 
taken at the same time, but different wavelengths, were then visually compared to 
ensure their features lined up. 

5.5 Dopplergrams and Intensities 

Velocity in in /s and line center intensity in arbitrary units were calculated at each 
pixel using equations 2.2 and 2.3. Because of the line centering problem described in 
section 4.1, the line center velocity and intensity were calculated using the -125, —40, 
and +40 mA filtergrams. The actual offsets of these measurements from the average 
line center was about —85,0, and +80 mA. Each dopplergram and intensity image 
was given a time mark averaged from the time marks of the three filtergrams. 

Filter tuning errors and seeing, which deteriorated toward the end of the run, 
produced some poor quality filtergrams in the data set. This in turned caused some 
nonsensical velocity and intensity determinations. Pixel velocities that were greater 
in magnitude than 10 kin/s and intensities that were a factor of 10 greater than 
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average were replaced with (lie mean velocity/intensity for that frame. 

5.6 Image Drift 

The rotational regridding assumes the images were rotating about a central point. 
In fact they rotate about the point centered in the spot tracker. This inconsistency 
produces a translational motion of the solar image. The point (on the sun) centered 
in the spot tracker, moreover, changes with time. This is caused by variable weighting 
of parts of the active region as the region evolves and as seeing conditions change. 

From equation 2.27. il is evident that to keep oscillation frequencies from shifting 
by more than 60 //Hz (0.5 frequency bins) at a wavenumber of 3 Mm” 1 , the image 
drift must be less than 120 m/s. This is equivalent to 1.4 pixel/hr. Image drift is 
measured from the cross-c orrelation of successive images. (It is defined as the distance 
of the maximum of the cross-correlation from the origin.) 

The image drift is traced out over time in figure 5.2. It was measured using the 
central 64 x 64 pixels of every tenth —125 mA filtergram. As can be seen from the 
figure, drift can be measured to significantly better than 1 pixel. The measured drift 
was smoothed by convolving the time series of each component of the vector with a 
Gaussian of 6 bin fwhm. The trend was then interpolated to the dopplergram and 
intensity time marks and removed from these images by regridding. 

5.7 Angular Scale Determination 

For full disk images of the sun, the image scale can be determined by measuring 
the radius of the solar disk. For high resolution images, however, the full disk is 
not visible, and indirect methods must be used. The scale of the La Palma images 
(0.8x lens) was determined by comparison of magnetograms from La Palma and the 
Kitt Peak National Observatory (figures 4.4 and 4.2). The Kitt Peak magnetogram 
was made at 15:17 UT on August 2, 1990, and has a scale of 1 "/pixel to within 1% 
(Harvey 1991). The La Palma magnetogram was interpolated from magnetograms 
made at 13:35 UT and 16:27 UT. 
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Figure 5.2: Image drift during observing run. Image motion determined by cross- 
correlation is traced out for the observing run. Points represent filtergrams taken 177 
s apart; solid line shows low-pass filter of points; circles represent 18 min time marks 
(< = 0 is at upper right). 




CHAPTER 5. DATA REDUCTION 


70 


The interpolation method accounted for evolution of the active region by measur- 
ing the motion of the magnetic field structures using local correlation tracking. The 
two magnetograms were then stretched by the time-weighted local displacement grid. 
Finally, the time- weighted signals of the stretched magnetograms were added. 

To compare magnetograms, 43 points on the La Palma magnetogram were flagged 
manually by cursor on a PIXAR computer. Corresponding points were flagged on the 
Kit t Peak magnetogram. Distances between a point and the centroid of the 43 points 
were then compared for the two magnetograms. The field of view for the La Palma 
images was, thus, determined to be (224.39 ±0.70) Kitt Peak pixels. Assuming the 
Kitt Peak scale measurement . 

La Palma Field of View = 224” ±2” (square) (5.3) 
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Chapter 6 

Measurements from Power 
Spectra 


The sets of dopplergrams and intensities were apodized in the space and time di- 
rections with Blackman- Harris functions (equation 2.32). Three-dimensional FFT’s 
were calculated for the full field of view, as well as for the three 79" x 79" subfields 
shown in figure 6.1. The temporal and spatial Nyquist frequencies were 28.3 mHz and 
9.9 Mm -1 , respectively. The frequency bin size was 117 /iHz, and the wavenumber 
bin size, 0.039 Mm -1 for the full field, and, 0.11 Mm' 1 for the subfields. 

The squared modulus of each Fourier transform was calculated to produce a 3-d 
power spectrum P(k',v). The dispersion relations given by equations 1.7 and 1.9 
are symmetric about the k — 0 axis. They define “trumpet” shaped surfaces (Hill 
1988) of approximate form v oc |£| 1/2 . (See figure 6.2.) This symmetry allows us 
to integrate the power over the azimuthal coordinate, in order to increase the signal 
to noise ratio. The resultant, collapsed power spectrum P(k,is) will show ridges of 
power at frequencies given by the dispersion relations for the p-modes and /-modes. 

Properties of the oscillations, such as the dispersion relation and mode power, were 
measured from the collapsed power spectra. This was accomplished by fitting one- 
dimensional slices of the spectra. Fits were then compared with theoretical predictions 
to evaluate models of the solar atmosphere. Magnetic versus nonmagnetic regions 
were compared to search for differences in wave properties between the two regions. 
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Figure 6.1: Magnetogram of AR61SL The three outlined subfields were used sepa- 
rately to measure wave properties in network field and quiet regions. The small tick 
marks are 5". 
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6.1 Collapsing 3-d Power Spectra 

The azimuthal symmetry of equations 1.7 and 1.9 is realized only in solar coordinates, 
comoving with the gas. In practice, one or both of these conditions may not be met. 
The frequency shift described by equation 2.27 occurs, with a cosine dependence on 
azimuth 9 , if there is a relative (horizontal) motion between the observer’s coordinates 
and the gas. The rigid alignment procedure described in section 5.6 keeps the observa- 
tional coordinates fixed with respect to the sunspot. It has been shown, however, that 
there is relative motion between magnetic flux and nonmagnetized gas (Snodgrass Sz 
Ulrich 1990). Hill (1988) finds (lows of ~ 100 m/s below the photosphere. 

Visual inspection of the power spectra did not reveal evidence of flows. Thus, 
if present, flows produce frequency shifts of less than one frequency bin (117 //Hz). 
If plane waves propagated isotropically, a sinusoidal frequency shift with 9 would 
serve to widen the ridges in the collapsed power spectra. In the measured spectra, 
however, power depends mildly on angle 0. Thus, a horizontal flow would produce a 
slight asymmetry in addition to broadening a ridge. 

The second condition affecting the symmetry of the trumpets is geometrical. For 
observations away from disk-center, the projection of solar coordinates onto obser- 
vational coordinates causes foreshortening. Waves propagating toward or away from 
disk-center appear with shorter wavelength (higher wavenumber) than waves prop- 
agating in the perpendicular direction. This foreshortening distorts the cylindrical 
symmetry ol the trumpets to give them an elliptical rather than circular cross section. 
Procedures to measure and compensate for this effect are described below. 

6.1.1 Distance Scale Determination 

The distance scale of the images or, equivalently, the wavenumber scale of the power 
spectra follows from trigonometry, once the magnitude and direction of foreshortening 
have been determined. This measurement was made by fitting an ellipse to the cross 
section of a. trumpet in a 3-d power spectrum. 

A constant frequency slice, in the 5-minute band, of the full-field Doppler spectrum 
is shown in figure 6.3. The slice has been remapped to cylindrical polar coordinates. 
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Figure 6.3: Constant frequency slice ot the 3 - cl power spectrum. A slice at 3.99 mHz 
has been remapped to cylindrical polar coordinates. Log [power) is shown as a function 
of wavenumber and azimuth. 

So, the cross section ot a trumpet appears at roughly constant wavenumber for all 
values of 0. These horizontal bands of power are the "rings" described by Hill (1988). 
In polar coordinates, an ellipse centered at the origin has the form, 

It = A- 0 ['l - (1 - /x 2 ) cos 2 (0 -0 o )]~i (6.1) 

where k 0 is half the minor axis (the true wavenumber); 0 Q is the angle of the major 
axis (direction toward disk-center); and p is the ratio of minor to major axes (the 
fractional foreshortening). A close 1 inspection of figure 6.3 shows that each of the 
bands exhibits this oscillatory behavior ot k as a function of 9. 

The values of k 0 .() 0 . and ft were fitted to the measured wavenumber k{0) for the 
Pi ring in the frequency slice 1 u = 3.99 mHz. The function k{0) was defined as the 
wavenumber at peak power in the constant frequency slice. It was determined by 
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quadratic interpolation. Ihe measured fractional foreshortening was, 

/< = 0.956 

This corresponds to a distance from disk-center of 17°, which agrees with an approx- 
imate measurement made from the full-disk magnetogram in figure 4.2. 

The assumption of a planar atmosphere has been made here, in treating the fore- 
shortening as constant across the field of view. If the foreshortening varies appreciably, 
modes of different wavenumber will be mixed in the 3-d power spectrum. This will 
smear the trumpet surfaces in (A \v) space in a manner which cannot be corrected. 
For large fields of view, the images must be mapped to "flattened” coordinates before 
transforming to Fourier space. The aplanarity for the data analyzed here was not 
great enough to warrant such treatment. 

The 224" field of view represents 13° along a great circle on the Sun. Since the 
average distance from disk-center was about 17°, plane waves propagating along the 
direction of foreshortening will be mixed over a wavenumber range of ±3.5%. While 
this may seem large, only the waves propagating exactly along a diameter of the disk 
experience mixing to this extent. Waves propagating in the perpendicular direction 
are not mixed at all. The net widening of the ridges is, thus, small. This can be seen 
in figure 6.3. The width of the ridges appears independent of 0. The width, therefore, 
is not an artifact of the coordinate system. 


6.1.2 Azimuthal Projection 

The power was integrated around ellipses in the {k x ,k y ) plane. To do this, A' 0 was 
computed from equation 6.1 for each ordered pair {k x ,k y ) in the 3-d spectrum, and 
power from that bin was added to the appropriate k bin in the 2-d spectrum. Bins 
in the 2-d spectrum were normalized by (multiplying by) the ratio of the area of an 
elliptical band to the number of bins accumulated from the 3-d spectrum. The area 
A of an elliptical band was approximated by, 

A — tc— \/ci 2 + b 2 


( 6 . 2 ) 
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Figure 6.1: Scaled velocity- power spectrum. A portion ot the collapsed power spec- 
trum for the full field of view, with variable gray-scaling in the k and v directions. 
Background has been subtracted. The noise visible at the top of the spectrum is an 
artifact of the scaling. From Fernandes 1 1 al . (1992). 

where a and b are the major and minor axes ot the ellipse; and w is the width of a 
wavenumber bin in the 2-d spectrum. 1 he two-dimensional spectrum P(k\i/) for the 
full-field dopplergrams is shown in figure 6.-1. 
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6.2 Parametrization and Least-Squares Fitting 


6.2.1 Fitting Procedure 


The frequencies, widths, and heights of the ridges were measured as functions of 
wavenumber. To do this, constant wavenumber slices of the power spectra were indi- 
vidually fitted by least-squares method. Examples of slices at constant wavenumber 
appear in figure 6.5. A single isolated mode should have a Lorentzian profile. 


/» 


PqF /2tt 

(p - vq) 2 + r 2 /4 


(6.3) 


where P 0 is the mode power: / is the mode frequency; and T is the fwhm. The high 
frequency mock modes described by Kumar et ai (1990) should have a sinusoidal 
profile. 

sm 2 (a,\r 0 /c) 


P{v) -x 


(6.4) 


jj 2 — u; c 2 sin 2 (u;a/c) 

where .r 0 is the source depth; jj c is the acoustic cutoff; and a is the height of the upper 
reflecting surface. Solar processes such as horizontal flow fields, and measurement 
effects such as seeing and apodization tend to smear the shape to produce a Gaussian 
profile. 


P(i') 


Po 




(6.5) 


<7\/Z7r 

where P 0 is the mode power; is the mode frequency; and a ~ fwhm/ 2.36. 

The ridges lie on top of a background which decreases monotonically with wave- 
number. In the photosphere, the granulation is caused by overshoot of convective cells 
into the overlying convec lively stable layers. The power, in the photosphere, below 
the /-modes is caused by this convective overshoot. Extrapolation of calculations of 
Nesis N: Mattig (1989) to t he low chromosphere suggests a similar explanation for the 
background power there. 

Because of the uncertainty of the physical description for the background, there 
is no proper way to parametrize it. Nevertheless, the background must be accounted 
for to avoid contaminating t he ridge fits. It can either be included as an extra term in 
the fits, or fit independently and subtracted from the power spectrum before fitting 
the ridges. 
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Figure 6.5: Constant wavenumber slices of the 3-d power spectrum. Data are plotted 
as points with error bars. Background and individual Gaussian modes are plotted 
as dashed lines; and the sum of background and modes is plotted as a solid line. 
Top left: k = 1.16 Mm -1 ; Bottom left: k = 5.25 Mm -1 ; Right (top k bottom): 
k = 2.34 Mm' 1 . 





CHAPTER 6. MEASUREMENTS FROM POWER SPECTRA 


80 


The latter method was chosen because the background could be well characterized 
using the range of frequencies on eit her side of the ridges. The shape of the background 
Pk suggested a power law parametrization. 

Pk{y) = A k ,r B * ( 6 . 6 ) 

where Ak and B are parameters fitted independently for each wavenumber slice. 
The frequency ranges used for the background fits include: (0.12 < yj mHz < 

1.3) and (12.0 < y/xwWz < 15.0). 

The ridges were fitted after subtracting the background. As can be seen from 
figure 6.5, the ridges overlap, and so must be fitted simultaneously. With three 
parameters per Gaussian, and up to eight visible ridges, the fits involved up to 24 pa- 
rameters at a time. This large parameter space combined with noisy data necessitated 
the development of a very robust fitting algorithm. 

The fitting program is described in the next section. It provides a nonlinear, 
weighted, least-squares fil to ungridded data. The weights used in the least-squares 
fit were the measurement errors. Hence, the fitting procedure minimized \ 2 . The 
error on the power was estimated from equation 2.38. This in turn assumed the 
velocity and intensity errors derived in section 2.1. 


6,2.2 Levenberg-Marquardt Method for Nonlinear Regres- 
sion 


The basic algorithm follows the Levenberg-Marquardt procedure (Press et al. 1989). 
The method attempts to minimize the \ 2 in parameter space (a), given N measured 
values yi ± a t at points .r,. 





cr t 


2 


(6.7) 


where //(.»*: <7) is the Ulted function. 

The minimum of the surface \ 2 {a) is searched for in an iterative manner. Succes- 
sive guesses for the parameter vector a are made based on the gradient and curvature 


CHAPTER 6 . MEASUREMENTS FROM POWER SPECTRA 


81 


matrix for the \ 2 surface. When the distance to the minimum, as estimated from 
the curvature, is large, the gradient is weighted more heavily in the guess. When the 
minimum is near, the curvature carries more weight. The algorithm assumes it has 
found the minimum when it can no longer reduce \ 2 by more than 0.1. The covariance 
matrix for the fit is simply the inverse of the curvature matrix at the minimum. 

To prevent runaways, the parameters are tested after each iteration of the algo- 
rithm in which \ 2 is reduced. Parameters that are not within set limits are fixed at 
appropriate values, and no longer varied in the fit. The remaining parameters are 
reset to their initial value's o {) . and the procedure starts again from the beginning. 

For Gaussian ridges, the limits for the height parameter are e (small, positive) 
and h inax (very large). Value's falling outside the limits are fixed at the limit. The 
limits for the width parameter are one-half and twice the expected width. Again, 
values falling outside the limits are fixed at the limit. The limits tor the frequency 
parameter fall hallway between the expected positions for ridges ( no ) • For the bottom 
(top) ridge in a fit, the lower (upper) limit is defined such that the two limits are 
equidistant from the expected ridge position. Values falling outside the limits are 
fixed at the initial guess ff 0 . 

The reason for this treatment of parameter limits is to prevent noisy ridges from 
disturbing the fits of adjacent ridges. Ridges which are too noisy to be properly fitted 
are treated as background . 

For the first spectrum to be fitted, initial guesses at a number of wavenumbers were 
made by eve, using cross hairs on a workstation display. Guesses at other wavenum- 
bers were interpolated or extrapolated from the nearest four fitted wavenumber slices. 
For the remaining power spectra, initial guesses were interpolated from the already 
fitted spectra. 

In addition to the background and ridge fits, the least-squares fitting program was 
used to interpolate and extrapolate mode frequencies to make the initial guesses, and 
to determine the parameters of equation 6.1 in section 6.1.1. 
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Figure 6.6: Power law exponent for background fit. Exponent vs. wavenumber. 

6.3 Results of Fits 

6.3.1 Background 

The exponent to the power law fit is plotted in figure 6.6. The background velocity 
power, integrated over frequency, is plotted in figure 6.7. The spectrum falls off much 
more slowly with wavenumber than the spectrum of convection in the photosphere re- 
ported by Chou et ai (1991 ). The root- mean-square (RMS) velocity is about 820 m/s. 
Nesis Mattig (1989) have measured RMS velocities, due to “mechanical overshoot”, 
that increase with height in the upper photosphere. Their model predicts the height 
dependence, for line-of-sight velocities, shown in figure 6.8. A smooth extrapolation 
of this curve to the height of formation of the Na Di line (500-700 km) yields an 
RMS velocity of 700-1100 m/s, in agreement with our measurements. 
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6.3.2 Mode Frequencies 

Frequencies ol ridges from the lull-field velocity and intensity spectra, and the quiet 
subfield velocity spectrum are tabulated with their respective statistical errors in 
appendix A. The velocity points are plotted in figure 6.9. All tabulated points 
correspond to fits in which the offset, width, and height of the Gaussian converged 
to physically reasonable values, and the integrated power under the Gaussian was at 
least 4 standard deviations above the background. 

All six fitted p-mode ridges extend up to at least 7.0 mHz. The ps and p$ ridges 
extend up to a frequency of 9.5 mHz. There is no deviation from a smooth extrap- 
olation of the ridges at lower frequencies. Libbrecht, Woodard Sc Kaufman (1990) 
have made measurements that overlap with ours at frequencies up to 5.3 mHz. As 
can be seen from figure 6.10, the fitted frequencies are consistently higher than those 
of Libbrecht et al. by an average of (35 ± 13) /xHz. This difference, however, is only 
a quarter of the frequency bin width, and is indicative of the systematic error in the 
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Figure 6.8: Calculated Doppler velocity vs. height for the Fe I (5123.75 A) line. From 
Nesis L Mattig (1989). 

measurement. 

The /-ridge is visible at wavenumbers up to 5.57 Mm -1 . There is, however, a 
systematic discrepancy between frequencies fitted from the velocity and intensity 
spectra as can be seen from figure 6.11. Frequencies measured from the velocity 
spectra are consistently higher than those from the intensity spectra by up to 50 /zHz 
at wavenumbers below 2.5 Mm -1 . This difference increases to ~ 100 /zHz at 4 Mm -1 . 
The discrepancy may in part be due to a difference between background shapes in 
the velocity and intensity spectra. 

Libbrecht et al. (1990) report a Gaussian background, centered near 3 mHz, in 
their velocity spectrum. Inclusion of such a background in our velocity fit would bring 
the /-ridge frequencies into better agreement with the intensity data at wavenumbers 
below 1.1 Mm -1 , but would push them farther apart at higher wavenumbers. Fitting 
such a background function, moreover, would require a longer time series with corre- 
spondingly better signal-to-noise ratio. It is noteworthy that the /-mode frequencies 
measured from the intensity spectra differ with those reported by Libbrecht et al. by 
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Figure 6.9: Measured mode frequencies (from velocity spectra): v vs. log 10 k. Fre- 
quencies of the / and p,-p (i ridges are plotted with 1 a error bars. For k < 3.82 Mm' 1 , 
frequencies are taken from the full-field spectrum. At higher k, values from the quiet 
subfield are plotted. 




86 


CHAPTER 6. MEASUREMENTS FROM POWER SPECTRA 



0.4 0.8 1.2 1.6 2.0 2.4 

k (mM* 1 ) 

Figure 6.10: Comparison of measured p-mode frequencies with Libbrecht et al. (1990). 
Frequency difference vs. k is plotted for Pi,P 2 ,P 3 (from top down). 
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Figure 6.11: Comparison of /-mode frequencies measured from velocity and intensity. 
Frequency difference vs. A*. 

under 30 ^iHz, as shown in figure 6.12. 

The weighting of magnetic regions is different for intensity and velocity spectra. 
However, as shown in section 6.3.5, the discrepancy is not caused by a difference in 
mode frequencies in the magnetic and nonmagnetic regions within the field of view. 
Regardless of which set of /-mode frequencies we use, the measured frequencies lie 
below the theoretical frequencies given by equation 1.9 at wavenumbers greater than 
1 Mm' 1 . The difference is 100-200 /iHz at 4 Mm' 1 . The fractional difference between 
measured and theoretical frequencies is plotted in figure 6.13. 


6.3.3 Mode Widths 

Ridge widths are plotted as error bars to the frequencies in figure 6.14. The ridge 
widths at low frequency are due to the apodizing function. At high frequency, the 
ridges are broader than expected for mock modes. The measured widths probably do 
not represent the true widths. The width is the most sensitive parameter to noise in 
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Figure 6.12: Comparison of /-mode frequencies from full-field intensity spectrum with 
Libbrecht et al. (1990). 

the fitted power spectrum. 

6.3.4 Mode Power 

Mode power is defined as the integral under the fitted Gaussian. 

p = (6.8) 

where h is the height of the ridge; and w , the fwhm. The error on the power is given 

by, 



where cr^w 2 is the covariance of the height and width. The mode power is plotted 
as a function of frequency in figure 6.15. The p-mode power has the same frequency 
dependence for all p-ridges. There is very little explicit wavenumber dependence. 
Power increases as u 73 to a maximum at 3.4 mHz, at which point it turns over and 
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Figure 6.13: Comparison of measured /-mode frequencies with theory. Frequencies 
from the full field intensity spectrum are compared with the relation u> 2 = gk (g = 
274.1 m/s 2 ). 
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Figure 6.14: Frequency and ridge width vs. k. Top plot shows /-modes and even 
order p-modes. Bottom plot shows odd ordered p-modes. Error bars represent the 
fwhm from the fit. 
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falls off again as i/~ 4 6 . The /-modes have less power at a given frequency than the 
p-modes. As surface grav ity waves with very little compression they are physically 
different from the p-modes and are not expected to show the same power spectrum. 

An MTF correction was made to the power plotted in figures 6.7 and 6.15 by com- 
bining the calculated MTF of the telescope and CCD with the estimated atmospheric 
MTF. See sections 3.4 and 5.4.1. 

6.3.5 Comparison of Magnetic and Nonmagnetic Regions 

Mode frequencies were compared for the quiet, canopy, and sunspot subareas shown 
in figure 6.1. The frequency differences are plotted in figures 6.16 and 6.17. The 
frequencies of ridges in each subarea agree with those of the other subareas to within 
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Figure 6.16: Mode frequency differences for canopy and quiet subfields. Top left, 
/-mode; top right, p x ; bottom left, p 2 ; bottom right, p 3 . 

the scatter of the data (±40 /tHz). 

The experimental scatter is much larger than the error bars derived from the fit. 
This may simply indicate an underestimation of the velocity error, or it may be a 
result of solar noise. Modes of finite lifetime have an intrinsic frequency uncertainty. 
Determination of the frequency of the cavity resonance i/o requires observations for 
several mode lifetimes. The 2.4 hr observing run probably constitutes less than one 
mode lifetime over much of the wavenumber range. Hence, the measured frequencies 
will be distributed about uq in a manner which depends on the mode lifetime, and 
shape of the wavepackets. 
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Figure 6.17: Mode frequency differences for sunspot and quiet subfields. Top left, 
/-mode; top right, pi; bottom left, p 2 ; bottom right, 






Chapter 7 


Spatial Distribution of Power 


The spatial distribution ol power was measured tor groups ot modes, as well as for the 
sub- fundamental region ot the power spectrum. I'o do this, the set of dopplergrams 
V’(x, t) was masked using 11 (.r) trom equation 2.33, and a 3-d FFT, calculated for the 
full field of view. A spectral mask .\/(A\//) was then applied to the Fourier transform 
to remove all but a small range of wavenumbers ot a given ridge. An inverse Fourier 
transform of this masked 3-d spectrum yielded velocity 1 t ) for a specific set of 
modes. 

\\,(.?.t) = ^- l {.\I(h.^)^{W(x) V(.r,0}} (7.1) 

The mean power was calc ulated at each pixel to make a velocity power map P\f(x). 

P.m(. r) = 2 / r |V A/ (.r.O|" (It (7-2) 

I Jo 

= 2 f' vq |^-'{-U(C/,)^ 3 {U'(.r) mo}}| 2 dv (7.3) 

where JT 7 1 denotes the spatial, inverse Fourier transform. 

7.1 Mode Masks 

Masks were constructed fur the / through p :r modes in 6 wavenumber ranges. In addi- 
tion, two sub-fundamental and a high trequencv set of masks were made for the same 
wavenumber ranges. Tin* mask limits are shown in figure 7.1. Three-dimensional 
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masks were made from these 2-d limits by rotating the loci in the azimuthal direc- 
tion. To account for the foreshortening, the rotation followed ellipses described by 
equation 6.1. 

7.2 Results 

Two sample power maps are compared with a magnetogram and an H a image in 
figure 7.2. Only the central disks of the frames are visible due to the apodizing 
window. The sub^/and /-mode maps were both made for the wavenumber range 
1.95 < kj Mm -1 < 2.31. Ibis is equivalent, to a wavelength of about 4", or 3 Mm. 
The mottled appearance of the power maps is an artifact of random interference of 
waves in a finite time series. 'The sub i-/ map shows little variation other than an 
overall trend of higher power in the upper right hand corner. The /-mode map shows 
a depression in power in the active region. 

A selection of power maps is shown in figure 7.3. An estimate of the RMS velocity 
in each map can be obtained from figure 7.4. The power in high wavenumber modes 
is depressed in the active* region. At low wavenumbers, the active region has little 
effect on the power. The snip -/and sub 2 -/maps show no correlation of power with 
magnetic field at any wavenumber. The high frequency maps, though not shown in 
figure 7.3, looked very similar to p-mode maps of the same wavenumber range at 
lower frequency. There is little difference between the /, p 1? p 2 » and P 3 maps for a 
given wavenumber range. 

A horizontal slice through the /-mode map from figure 7.2 was made above the 
sunspot. The slice was low pass filtered using a Gaussian impulse response function 
(16 bin fwhm). See figure 7.5. This low pass filtered slice is compared with the 
magnetogram signal in figure 7.6. The power depression is significantly broader than 
the spatial extent of the photospheric magnetic field. It is narrower than the canopy 
field at ~ 1200 km, as can be seen from the H a image in figure 7.2. The filaments 
visible in II a are streams of gas traveling along magnetic field lines. The filaments, 
thus, trace out the field lines which form a canopy above the photosphere and lower 
chromosphere. 
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Figure 7.1: Masks for power distribution maps. Measured mode frequencies are 
shown as dots. Mask limits are shown as solid lines. The high frequency (HF) masks 
covered the range 6.0 < vf mHz < 10.0; the p-mode mask limits were set halfway 
between modes; the /-mode limits were set equidistant from the /-mode; the sub 2 -/ 
band extends from the /-mode down to u = 0; the sub x -/band is the same width as 
the /-mode band. It is a subset ot the sub 2 -/band. 
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Fi gun' 7.2: Comparison ol power maps with magnetic fields. Top left: log power 
for sub | -/ ( 1 .9.' > < /.•/ Mm 1 < 2.51). lop right: absolute value of magnetogram (log 
scale). Bottom left: 11 n — IU0 mA image. Bottom right: log power for /-mode 
(1.95 < A/Mm -1 < 2.51). 1 he small lick marks are 5T 
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Figure 7.3: Power (list rilmt ion for selected masks. Log power is displayed for sub 2 -/. 
f, p i, pi- and />;j maps with wavenumber ranges: 

[1.95.2.31], [1.17.1.50]. [11.78.1.17]. [0.51 . 0.78] Mm" 1 . 
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Figure 7.4: RMS velocity vs. k. Curves fitted from quiet subfield velocity spectrum. 
Curve labeled sub-/ is taken from the background fit. 

It is surprising that t lie sub-fundamental power is not correlated with magnetic 
field. The RMS convective power in the photosphere is a factor of 2 lower in plage 
(network field) than in nearby quiet regions (Title et al. 1992). Moreover, the sub 2 -/ 
band lies just below the /-mode band, but shows completely different behavior. At 
the same wavenumber, and almost the same frequency, oscillation amplitudes are 
reduced in active regions, whereas non-oscillatory motion is not. 

It is impossible to distinguish between instrumental backgrounds and solar back- 
grounds from a single data set. It is difficult, however, to come up with an obser- 
vational effect that would produce the measured background. Filter transmission 
variation should affect only the very low wavenumber part of the spectrum. Seeing 
variations should be independent of frequency; whereas, the observed background 
decreases with frequency. 
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Figure 7.5: Slice of /-mode map. Horizontal slice from the [1.95,2.34] Mm 1 map. 
Also shown is a low pass filter of the slice. 
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Figure 7.6: Comparison of /-mode slice with magnetic field. Dashed line shows low 
pass filter of slice (arbitrary scale). Solid line shows magnetogram signal (arbitrary 
scale). 



Chapter 8 

Conclusions and Discussion 


8.1 Surfaces and Cavities 

The existence of a chromospheric cavity is expected to produce a resonance consist- 
ing of a horizontal band of power and "avoided crossings” at about 5.5 mHz. (See 
figure 1.4.) This cavity is required by the model of Balmforth & Gough (1990) to 
explain the existence ol high frequency modes. The /;-mode and /-mode ridges shown 
in figures 6.4 and 6.9 arc* clearly visible at 5.5 mHz, but there is not a horizontal band 
of power nor are there "avoided crossings". The clear absence of this resonance in our 
data shows that there is no chromospheric cavity. Either the reflection coefficient at 
the base of the corona is loo low or the surface is too irregular to form a good cavity. 

This result supports the alternate hypothesis of Kumar et al. (1990), that the 
high frequency ridges are a purely geometrical artifact. The result is in agreement 
with the conclusion of Woodard Libbrecht (1991); and it contradicts the claim of 
Dame et al. (1984) that there is a horizontal ridge of power at 5.5 mHz. Our power 
spectrum has a significantly better signal to noise ratio than that of Dame et al. 

The /-modes are sensitive to conditions in the chromosphere and transition re- 
gion. The fact that /-mode frequencies are lowered at high wavenumber may indicate 
interaction with chromospheric magnetic fields, energy leakage into the corona, or 
horizontal inhomogeneities in the chromosphere- corona boundary. The last of these 
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three possibilities is interesting from the standpoint of acoustic cavities. Perhaps ir- 
regularities in the interface between the chromosphere and corona are responsible for 
both lowering the /-mode frequencies, and for reducing coherent reflection of acoustic 
waves. 


8.2 Waves in Magnetic Fields 

There is no frequency shift for / and p-modes in magnetic network fields of 100 — 
— 1000 gauss to within the scatter ot t lie data (~ 50 //Hz). This limit is not stringent 
enough to determine whether the frequency shifts measured by Woodard et al. (1991) 
residt from small shifts over the whole disk, or larger local shifts, in active regions, 
which contaminate the full-disk measurement. The modes measured by Woodard et 
al. have widths of order 1 //IIz. Hence, to contaminate such a measurement, mode 
frequencies in active regions would have to be distributed within a few //Hz of the 
quiet sun values. A measurement of such precision would require observations of at 
least 20 hr duration, assuming the waves are coherent for the observing period. At 
La Pal ma, 20 hours of observat ions can be made in two days. Hence, the waves would 
add coherently if the mode lifetime is > 34 hr. If the mode lifetime is less than this, 
the signal to noise ratio increases only as the square root of time, so a much longer 
observing time would be necessary. 

Power maps for oscillations in the low chromosphere show that, at high wavenum- 
ber, wave amplitudes arc' suppressed 15 Mm outside of the photospheric magnetic 
field. (See figure 7.6.) In the photosphere, the power is suppressed within one (hori- 
zontal) wavelength from the field (Tarbell et al. 19SS). The vertical distance between 
the photospheric and chromospheric observations is under 0.5 Mm, and the (true) 
wavelength of the oscillations is ~ 2 Mm. Diffraction cannot account for the spread- 
ing of the region of the power depression because the waves have traveled less than 
one wavelength in the interval over which the spreading occurs. The power suppres- 
sion must, therefore, be produced by interaction of the oscillations with the field in 
the layer in where the warts are observed. 

The area containing magnetic flux was not measured directly at the height of the 
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observations. Its extent can be inferred, however, from the magnetogram made in 
the mid-photosphere, and the II o image in the mid-chromosphere. The field lines 
spread out rapidly above the photosphere. Hence, the extent of the field at the 
observed height should bo intermediate between the limits of the mid-photospheric 
and mid-chromospheric holds. This is consistent with the above conclusion. 

The /-modes propagate without compression, so they produce no pressure fluctu- 
ations in the region where oscillations are suppressed. The mechanism suppressing 
the wave motion must, therefore, act through the velocity fluctuations. 

8.3 Directions for Future Work 

The frequency measurements presented here show the existence of high frequency and 
high wavenumber oscillations. They show general characteristics of the dispersion 
relations such as the smooth continuation above the acoustic cutoff. They are not, 
however, high precision measurements in comparison with existing low wavenumber 
observations. Precision ol the' measurements can be improved with longer observing 
times. This would improve the signal to noise ratio, and enable better fitting of the 
background and ridges. \\ ith better frequency resolution, it would become necessary 
to correct the fitted frequencies for their weighting with mode power (Korzennik 
1991). 

Important characteristics ol the solar oscillations can be inferred from their phases 
(Deubner k Fleck 1989, 1990; Fleck k Deubner 1989). Trapped and nontrapped 
oscillations can be distinguished by the phase difference between velocity and intensity 
measurements. Propagating and evanescent waves can be distinguished by the phase 
lag between velocities at different. heights in the atmosphere. An attempt was made 
to measure velocity-intensity phase differences. The attempt failed, possibly because 
of the large height range included in the measurement (see section 1.4), and possibly 
because of the large number of modes included in a single frequency bin (Braun et 
ai 1992a). 

The interaction of magnetic fields with oscillations can be better quantified by 
measuring both wave amplitudes and magnetic fields at a range of heights in the 
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atmosphere. Improvements that, have been made to the tunable filter system since 
the 1900 observing season will allow Doppler measurements restricted to either the 
core or wings of the Na l) t line. In addition, magnetograms can be made in the core 
or wings of the line by finding the difference between dopplergrams made with left 
and right circularly polarized light. 

Braun et at. (1992b) report halos of enhanced high frequency acoustic power 
around active regions. A search was made for such halos in power maps made with 
the HF masks shown in figure 7.1. No enhancement of power was found. This does 
not contradict tin 1 result of Braun et aL. as they measured intensity fluctuations in 
the photosphere, in a lower wavenumber regime. It would be interesting to repeat 
the analysis with continuum brightness measurements for a field of view of at least 
6 arc minutes. 

The cause of the sub-fundamental background is, currently, unknown. Experi- 
ments should be done to discriminate between observational and solar backgrounds. 
For example, a set of liltergrams taken at fixed wavelength with a high cadence 
could be used to make artificial dopplergrams, and check the effect of seeing on back- 
ground. Another set of liltergrams taken at the nominal, observing wavelengths could 
be made, using integrated sunlight, to check the effect of filter-transmission variation 
on background. 

An attempt was made to find the oscillation sources described by Brown (1991). A 
temporally high-pass filtered movie of velocities was scanned visually for evidence of 
spatially and temporally localized sources of power. No good candidates were found. 
A better method might be to construct the two-dimensional acoustic Pointing vector 
described by Braun et at. (1992b). Large positive divergences of this vector field for 
high-pass filtered data could point out oscillation sources. 




Appendix A 
Frequency Tables 


Table A.l: Mode frequency v ± a u (/iHz) as a function of wavenumber k (Mm l ) — 
measured from dopplergrams. For k < 3.82 Mm 1 , frequencies are taken from the 
full-field spectrum. At higher k , values from the quiet-subfield are tabulated. 
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Table A. 2: /-mode frequency v ± cr u (/z Hz) as a function of wavenumber k (Mm l ) 
— measured from intensities. 


0.399 


0.437 

0.475 

1810 12 

1748 12 


0.513 


0.551 


0.589 


2089 11 


0.627 


0.665 

0.703 

0.741 

0.780 

0.818 

2176 11 

2228 11 

2281 11 

2318 11 

2373 11 


0.694 

0.932 

0.970 

1.008 

1.046 

2515 11 

2563 11 

2610 11 

2655 11 

2692 11 


1.064 


1.122 


2793 11 


1.160 


1.196 


1.236 


1.274 


2981 11 


1.312 

1.350 

1.388 

1.426 

1.464 

3023 11 

3066 11 

3103 11 

3147 11 

3173 11 


1.502 


1.578 

1.616 

1.654 

1.692 

1.730 

3302 11 

3336 11 

3361 11 

3412 11 

3438 11 


1.806 

1.844 

1.882 

1.920 

1.958 

3520 11 

3547 11 

3581 11 

3615 11 

3641 11 
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111 


1.996 


2.224 


2.034 


2.072 


2.186 


3740 11 

3773 11 

3811 11 

3842 11 


2.300 

2.339 

2.377 

2.415 

3931 11 

3967 11 

3996 11 



4026 11 


2.453 

2.491 

2.529 

2.567 

2.605 

4067 11 

4097 11 

4117 11 

4154 11 

4191 11 


2.643 


2.681 

2.719 

2.757 

2.795 

2.833 

4245 11 

4283 11 

4306 11 

4336 11 

4359 11 


2.909 

2.947 

2.965 

3.023 

3.061 

3.099 

4418 11 

4424 11 

4448 11 

4491 12 

4522 12 

4525 12 


3.137 

3.175 

3.213 

1 — 

3.251 

3.289 

3.327 

4535 12 

4564 12 

4596 12 

4646 12 

4668 12 

4660 13 


3.365 

3.403 

3.441 

3.479 

3.517 

3.555 

4716 13 

4726 14 

4812 13 

4832 14 

4828 14 

4902 14 


3.593 

3.631 

3.669 

3.707 

3.745 

3.783 

4874 15 

4910 17 

4914 16 

4948 16 

4896 20 

4961 20 


3.821 


5037 19 
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